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ABSTRACT 


This  thesis  presents  a  study  of  the  limiting  distributions  of 
response  probabilities  arising  from  two  stochastic  learning  models  derived 
by  Bush  and  Mosteller  [1]„  The  limiting  distributions  of  the  Experimenter- 
Controlled  Model  satisfy  a  functional  equation  which  has  proved  to  be 
extremely  valuable  in  determining  properties  of,  and  solutions  for  the 
distribution  functions.  Under  certain  restrictions  on  its  parameters, 
the  functional  equation  yields  a  direct  construction  of  purely  singular 
Cantor-like  functions.  Under  other  conditions,  absolutely  continuous 
functions  may  be  constructed  from  the  functional  equation.  This  leads 
to  a  sufficiency  condition  on  the  parameters  for  absolutely  continuous 
distribution  functions.  Necessary  and  sufficient  conditions  for  the 
limiting  distribution  to  be  the  uniform  distribution  have  also  been  found. 

The  limiting  distributions  of  the  Subject-Controlled  Model  satisfy 
a  Volterra-like  functional  integral  equation.  Due  to  the  complexity  of 
this  equation,  its  value  in  deriving  useful  results  has  been  rather  limited. 
However,  it  does  display  some  properties  of  the  limiting  distributions 
in  restricted  cases. 

The  last  half  of  the  thesis  is  concerned  with  methods  of 
calculating  approximations  to  the  limiting  distributions  by  numerical 
solutions  of  the  distribution  functional  equations.  An  error  analysis  of 
the  procedures  is  impossible  at  this  stage,  but  comparing  the  numerical 
solutions  with  known  results  in  special  cases  indicates  that  the  numerical 
procedure  is  sufficiently  accurate  for  most  practical  purposes. 
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CHAPTER  I 

INTRODUCTION 

Since  19^9  a  number  of  mathematicians  and  psychologists  have 
been  working  on  the  formulation  of  mathematical  models  to  describe  learning 
processes.  Thus  far,  these  models  have  been  restricted  to  the  rather  special 
situations  of  psychological  learning  experiments  because  of  the  simplifications 
imposed  by  experimental  controls. 

The  original  purpose  for  deriving  learning  models  was  to  present 
a  mathematical  framework  for  analyzing  statistical  data  from  psychological 
experiments.  Beside  being  of  practical  value  for  applied  statistics,  these 
models  have  given  rise  to  a  number  of  mathematically  interesting  problems. 

For  example,  some  new  problems  in  statistical  estimation  were  opened  up 
and  also  some  new  problems  in  the  solution  of  functional  equations  have 
arisen.  In  the  future,  learning  models  may  lead  to  new  insights  into  the 
processes  of  learning,  and  so  give  rise  to  new  and  better  learning  theories. 

The  basic  assumption  of  all  present  learning  models  is  that  the 
possible  responses  of  a  subject  can  be  considered  as  related  to  the  stimulus 
elements  of  the  environment  by  a  set  of  response  probabilities.  As  a  subject 
undergoes  a  number  of  experimental  trials,  its  responses  may  be  expected  to 
vary  as  learning  progresses.  This  change  of  behaviour  is  taken  as  a  reflection 
of  changing  response  probabilities.  Existing  learning  models  are  the  results 
of  attempts  to  describe  the  changes  of  response  probabilities  in  terms  of 
mathematical  transformations  and  to  deduce  the  consequences  of  these  trans¬ 


formations  . 
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One  of  the  earlier  models  is  described  by  Bush  and  Mosteller  in 
their  book  Stochastic  Models  for  Learning  [ 1 ] .  The  authors  begin  by 
defining  learning  as  "  ...  any  systematic  change  in  behaviour  ...  ", 
Furthermore,  learning  is  said  to  be  complete  when  some  type  of  stability 
of  behaviour  is  attained.  Using  this  operational  definition  of  learning 
gives  an  objective  basis  for  testing  the  model. 

The  general  Bush  and  Mosteller  learning  model  may  be  qualitatively 
described  in  terms  of  the  following  experiment  of  the  reward  training  of 
rats  in  a  T-maze.  A  T-shaped  maze  is  set  up  with  a  food  dish  at  each  of 
the  two  ends  of  the  cross-bar.  Each  experimental  trial  consists  of  setting 
a  hungry  rat  at  the  starting  point,  the  base  of  the  T,  and  allowing  it  to 
run  down  the  maze  to  the  choice  point  where  it  must  turn  either  to  the 
right  or  to  the  left.  Depending  upon  the  type  of  experiment,  the  rat  may 
find  food  or  an  empty  dish. 

The  mathematical  abstraction  of  the  learning  process  is  constructed 

from  two  considerations.  Suppose  that  at  the  n  th  trial  a  given  rat  has  a 

probability  p^  of  turning  right  and  a  probability  qn  =  1-Pn  of  turning 

left.  If  the  rat  turns  right  and  finds  food,  then  the  probability  of  it 

turning  right  on  the  succeeding  trial  will  be  increased.  Alternatively,  if 

the  rat  turns  right  and  is  unrewarded,  then  Pn+^  will  be  less  than  p^  . 

Bush  and  Mosteller  assume  p  ,  to  be  a  linear  function  of  p  such  that 

*n+l  *n 

U.l)  Pn+1  =  a  +  bpn  , 

where  a  and  b  are  parameters  depending  upon  the  individual  rat  and 


whether  the  rat  was  rewarded  or  not  on  that  trial. 
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From  this  general  model  of  varying  probabilities,  a  number  of 
special  models  have  been  constructed  to  represent  special  learning 
experiments.  The  problems  to  be  considered  in  this  thesis  arise  from  the 
"Experimenter-Controlled  Event  Model"  and  the  "Subject-Controlled  Event 
Model". 


The  Experimenter-Controlled  Model  may  be  illustrated  as  a  variant 
of  the  T-maze  experiment  in  which  the  experimenter  has  control  over  the 
frequencies  of  rewards  and  non-rewards  that  the  rat  recieves.  In  this 
learning  situation  the  experimenter  places  food  in  the  right  dish  on  a 
fixed  proportion,  ,  of  the  trials  on  which  the  rat  turns  right.  Similarly, 
he  places  food  in  the  left  dish  on  some  other  fixed  proportion,  >  of 

the  trials  on  which  the  rat  turns  left.  On  any  given  trial,  the  decision 
to  reward  or  not  to  reward  the  rat  is  made  by  some  suitable  random  procedure. 
The  rewarding  proportions  are  chosen  such  that  ^  =  1  . 

The  Subject-Controlled  Model  is  another  variant  of  the  general 
learning  model,  but  in  this  case  the  reward  or  non-reward  is  wholly  dependent 
upon  the  subject's  behaviour.  In  terms  of  the  T-maze,  the  experiment  is 
run  with  food  placed  on  the  same  side  at  each  trial.  If  the  rat  chooses 
the  reward  side,  it  receives  food,  otherwise  it  remains  hungry.  Thus  the 
outcome  of  each  trial  is  completely  subject  to  the  rat's  choice. 

Consider  the  effect  of  running  a  large  number  of  initially  identical 
rats  through  a  series  of  trials  of  a  T-maze  experiment.  Before  the  first 
trial  all  the  rats  possess  the  same  probability  p^  of  turning  right.  At 
the  end  of  the  first  trial  the  response  probabilities  of  each  rat  will  have 
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been  modified  depending  upon  the  direction  the  rat  chose  and  whether  or 
not  it  had  been  rewarded.  Thus  the  response  probabilities  of  the  rats 
will  now  be  distributed  amongst  several  values.  At  the  end  of  the  second 
trial  the  response  probabilities  will  again  have  been  modified.  The  value 
of  p^  will  depend  upon  the  value  of  p^  as  well  as  the  choice  and  out¬ 
come  of  trial  2.  Thus  for  each  value  of  p^  there  will  be  several 
possible  values  of  pg  .  As  the  number  of  trials  increases,  the  number 
of  possible  values  of  the  response  probability  increases  geometrically  so 
that  a  series  of  distributions  of  response  probabilities  is  generated. 

It  is  this  set  of  distributions  which  is  of  value  in  applied  statistics. 

Of  particular  interest  is  the  limiting  distribution  of  response  probabilities 
as  n  becomes  very  large  since  it  indicates  the  ultimate  behaviour  of 
the  rats  undergoing  the  learning  experiment. 

In  the  particular  cases  of  the  Experimenter-Controlled  and  Subject- 
Controlled  Models,  the  general  solutions  for  the  limiting  distribution 
functions  of  response  probabilities  are  unknown.  However,  it  will  be  shown 
that  these  functions  satisfy  certain  functional  equations  with  appropriate, 
auxiliary  conditions.  Unfortunately  there  are  no  known  methods  for  finding 
the  general  solutions  of  such  functional  equations,  but  these  equations 
are  useful  in  determining  some  of  the  properties  of  the  limiting  distributions, 
and  in  some  special  cases,  explicit  solutions  are  obtained. 

The  following  chapters  of  this  thesis  are  concerned  with  the 
properties  of  and  the  solutions  for  the  limiting  distributions  of  response 
probabilities.  In  chapters  II  and  III,  the  Experimenter-  and  Subject- 
Controlled  Models  respectively  are  derived.  From  these,  the  functional 
equations  which  are  satisfied  by  the  characteristic  functions  of  the  limiting 
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distributions  are  found.  The  characteristic  equations  are  then  inverted 
to  yield  the  functional  equations  satisfied  by  the  distribution  functions. 
The  remainder  of  Chapters  II  and  III  are  devoted  to  finding  some  of  the 
properties  and  solutions  of  the  functional  equations.  Chapters  IV  and 
V  are  concerned  with  the  problem  of  tabulating  numerical  approximations 
to  the  limiting  distributions  by  numerical  solutions  of  the  functional 
equations.  These  chapters  also  include  some  tables  of  the  numerical 
solutions  for  selected  cases.  The  final  chapter  summarizes  the  conclusions 


of  the  previous  chapters. 
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CHAPTER  II 

THE  EXPERIMENTER-CONTROLLED  MODEL 


§2.1  Introduction 

In  this  chapter  the  Experimenter-Controlled  Model  for  a  learning 
experiment  having  two  possible  responses  is  treated  in  detail.  The 
mathematical  formulation  of  the  basic  model  is  derived  using  the  assumption 
that  learning  can  be  considered  as  a  process  of  systematic  changes  of 
responses  to  given  stimuli.  Since  no  two  subjects  can  be  expected  to 
learn  at  exactly  the  same  rate,  there  are  many  possible  values  of  the 
response  probabilities  at  the  n  th  trial  of  the  experiment.  It  is  shown 
that  the  moment  generating  function  of  this  distribution  of  response 
probabilities  satisfies  a  certain  recurrence  relation.  As  the  number  of 
trials  becomes  large,  the  distribution  of  response  probabilities  tends  to 
a  limiting  distribution  whose  moment  generating  function  satisfies  a 
functional  equation. 

The  properties  of  this  limiting  distribution  will  be  the  major 
concern  of  this  chapter.  In  order  to  examine  this  distribution,  the  moment 
generating  functional  equation  is  transformed  to  a  functional  equation 
satisfied  by  the  characteristic  function.  A  further  transformation  is 
made  to  standardize  the  characteristic  functional  equation  before  It  is 
inverted  to  yield  a  functional  equation  satisfied  by  the  limiting 
distribution  function.  The  known  properties  of  the  limiting  distribution 
are  reviewed  and  a  number  of  additional  properties  and  explicit  solutions 
are  found  from  the  functional  equation. 
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$  2.2  Derivation  of  the  Model 

Consider  a  simple  learning  experiment  in  which  each  trial 
consists  of  presenting  a  subject  with  two  mutually  exclusive  alternatives 
and  A^.  The  subject's  choice  of  A^  is  called  response  ,  and 

the  choice  of  A^  is  called  response  .  Thus  R^  and  R^  must  be 

mutually  exclusive  and  exhaustive  responses. 

It  is  assumed  that  at  each  trial  the  subject  has  a  definite 
probability  of  choosing  each  of  the  two  alternatives.  These  probabilities 
are  denoted  Pr'Pi)  ,  (i  =  1,  2).  From  the  condition  of  mutual  exclusiveness 
if 

(2.1)  Pr^)  =  p  , 
on  a  particular  trial,  then 

(2.2)  Pr(R2)  =  1  -  p  , 
on  that  same  trial. 

In  order  that  the  experiment  be  meaningful,  it  must  be  arranged 
so  that  each  time  a  response  occurs,  an  outcome  follows;  let  0^  and  0^ 
be  exhaustive  and  exclusive  outcomes.  From  all  these  considerations  the 
following  events  can  be  defined; 


V  : 

Ri 

followed  by  0^  , 

S'  : 

R2 

followed  by  0^ 

(2.3) 

E2  : 

Ri 

followed  by  0 2  , 

F  •  • 

\  * 

R2 

followed  by  0 2 

Assuming  that 

the  events 

V 

and  E,  '  and  the 

4 

events 

V 

and  E  '  are 

3 

equivalent,  one  may  define  the  events 


■  . 

:}  ,  '  :  •  '*1  k<>30<-ab 
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(S.k) 


E1  :  V'-A'  • 

e2  :  E2'WIy  . 


Under  the  above  conditions,  the  new  events  and  are  exclusive 

and  exhaustive. 


The  basic  axiom  of  the  model  is  that  the  occurrence  of  E  (j  =  1,  2) 
on  any  trial  modifies  Pr'R^)  for  the  succeeding  trial.  In  the  interests 
of  simplicity,  a  linear  transformation  was  chosen.  Thus  if  Pr’R^)  =  p  on 


,  th  _  .  . 

the  n  trial 

,  then  the  value 

of 

Pr(R^)  on 

s  t 

the  n+1  trial  is  given 

by  applying  one 

of  the  linear  operators  or 

defined  by 

(2.5) 

V  -  ai 

+ 

(Xjp  , 

(2.6) 

V  =  as 

+ 

a2P  ; 

where  a^  and 

a.  (i  =  1,  2) 

are  parameters 

satisfying  the  conditions 

0  <  a.  < 

1 

9 

(2.7) 

—  i  — 

0  <  a.  < 

™  1  — 

1 

-  a. 

1 

The  choice  of  the  operator  is  determined  by  the  occurrence  of  event  or 

V 

In  experimenter-controlled  learning  experiments,  the  experimenter 
arranges  that  on  a  fixed  proportion,  x-^  ,  of  the  trials  on  which  R^  occurs, 
0^  follows.  On  the  remaining  trials  that  R^  occurs,  0^  follows.  From 
the  conditions  ^2. 3)  and  (2.^-),  this  means  that  x^  is  the  probability  that 
operator  is  applied,  and  x^  the  probability  that  is  applied. 


smsTsq  &?z  (f  i  *  i}  .,D  .  s  , 
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Thus 


(2.8) 


i*:^)  =  Kj  , 

Pr(E2)  =  *2  . 


If  the  subject  has  some  initial  response  probability  Pr *R^)  =  pQ  , 
then  the  possible  values  of  Pr(R^)  throughout  the  first  few  trials  will 
be  generated  as  shown  in  Figure  2.1. 


til 

Since  there  are  two  possible  events  at  each  trial,  at  the  n 
trial  there  must  be  2n  possible  values  of  Pr(R^).  Let  these  values  be 
denoted  by  p  (m=l,2,3> • • • >2n  J  n=l,2,3>»»»  )  •  Similarly,  let  the 

IQ  f  II 

corresponding  probabilities  of  the  occurrence  of  p^  ^  be 


(2.9) 


P 

m,n 


Pr  (p, 


m,n 


)  • 


Thus  if  Pr(R, )  =  p  on  trial  n  ,  then  on  trial  n+1  the  possible 

v  1'  rm,n 

response  probabilities  are: 


“l 

Pm,n 

ai 

+ 

alPm,n 

q„ 

P 

ao 

+ 

a„p 

2 

m,n 

2 

2  m,n 

with  probability  jt.  P 

r  J  1  m,n 


with  probability  P 

J  2  m,n 


(2.10) 
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TRIAL 


Figure  1 


0  1  2 \ 


Possible  ways  that  the  response  probablity  Pr(R^)  can  vary 
throughout  the  first  three  trials,  given  an  initial  value  pq. 
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The  upper  entry  of  each  pair  is  a  possible  value  of  Pr(R^) 
expressed  in  terms  of  the  operators  and  and  the 

initial  response  probability  pQ  .  The  lower  number  is  the 
probability  of  the  occurrence  of  that  particular  value  of  the 
response  probability.  (Reproduced  from  [1]  p.70). 

§2.5  Properties  of  the  Model 


For  the  purposes  of  statistical  analysis,  much  useful  information 

is  provided  by  the  moments  of  the  distributions  of  response  probabilities. 

til 

Although  there  is  no  known  general  expression  for  the  m  moment  at 
the  n^  trial,  some  recurrence  relations  for  the  moments  have  been  found, 
McGregor  and  Hui  [2]  derived  a  recurrence  relation  satisfied  by  the  moment 
generating  functions  by  using  the  following  argument. 


Let  mn(@)  be  the  moment  generating  function  of  the  distribution 


m  . (@)  we  have 
n+1 


m 


n+r 


at  the 

nfc^  trial.  By  the 

definition 

rH 

£  Loo’ll 

1! 

(p#  n+1  0) 

P,  e 

i,n+l 

9 

2n 

-T.  ■ 

i=i 

)0 

+  rt2Pi,ne 

aiQ 

=  V 

2n 

r-i  a.p.  6 

l  ♦ 

2n 
a20  V 

v  Lp 

(  ao+aoP  g  JP 


1=1 


1= 1 


Therefore 

(2.11) 


a„0 


a_0 


-j^  2 

mn+l(0)  =  *1®  mn(ai0)  +  *2e  mn^a20^  * 
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Obviously  the  moment  generating  function  of  the  distribution  of  response 
probabilities  at  the  n+1  trial  is  dependent  on  all  the  moment  generating 
functions  of  the  previous  trials. 

Of  particular  importance  to  psychologists  is  the  behaviour  of 
the  subject  when  complete  learning  has  been  attained.  Karlin  [3]  opened 
the  way  to  this  study  by  proving  that  the  limiting  distributions  exist 
as  the  number  of  trials  tends  to  infinity.  However,  relatively  little 
has  been  published  about  the  properties  of  these  limiting  distributions. 

It  is  this  problem  which  will  be  the  concern  of  the  remainder  of  this 
chapter. 

The  moment  generating  function  of  the  limiting  distribution 

of  response  probabilities  satisfies  the  limiting  form  of  equation  (2.11). 

Thus  if  lim  m  (@)  =  m(0)  ,  the  above  recurrence  relation  reduces  to  the 
n— ►  00 

functional  equation 

a  0  a  0 

(2.12)  m(0)  =  e  m(a^0)  +  e  ^(a^) 

We  transform  (2.12)  into  a  functional  equation  satisfied  by  the 
characteristic  function  flpp(0)  by  replacing  0  by  i0  giving 

ia  0  iao0 

(2.13)  Cpp(0)  =  e  +  *2  8  ^  'Pptag0)  * 

It  has  been  shown  by  Karlin  in  his  Trapping  Theorem  ([l]  p.98) 
that  the  limiting  response  probabilities  are  trapped  in  the  interval 
bounded  by  A^  and  A^  where 


(2.14) 


\  =  aiA1_ai) 


^2  —  * 


|  i.'B 


*  •;  fi  ■  '  ■ 
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Thus  if  F(p)  is  the  limiting  distribution  function  of  response  probabilities, 
the  Trapping  Theorem  implies  the  boundary  conditions 

F(p)  =0  ,  p  <  min  (A^  Ag)  ; 

(2.15) 

F(p)  =1  ,  P  >  max  (Aj,  Ag)  . 

Substituting  for  the  a^  (i  =  1,2)  in  equation  (2.15)  ,  we  obtain 

iA  (1-a  )0  iA  (1-a  )0 

(2.16)  9p(@)  =  ^e  cpp(a10)-Ht2e  ^  *  9p(a2@)  * 

For  the  sake  of  convenience,  a  change  of  variable  will  be  made 

to  map  the  limiting  variable  p  onto  the  interval  [-1,  1],  Without 

loss  of  generality  it  can  be  assumed  that  Ag  >  A^  since  interchanging 

A.  with  A  and  a.  with  a0  and  jt,  with  jt  in  equation  (2.16)  does 

1  d  1  d  1  d 

not  change  the  equation.  The  required  transformation  is 


(2.17) 


*  =  ((\  +  A2)  -  2p)/(Ax  -  Ag)  ,  (Ag  >  Al) 


By  definition,  9p(@)  =  E(e*^P)  ,  where  E(e*^P)  is  the  expectation  of 
e*^p  .  Making  the  transformation  (2.17)  then  yields 

i®[((^1  +  *2)  "  2P)/(A1  "  Ag)] 

cpx(0)  =  E(e  )  =  E(e 


=  e 


i@(A^  +  “  ^2^ 


9  (20/ (Ag  -  Ax)) 


Putting  x  =  20/(Ag  -  A  )  ,  we  get 


■is(A1  +  Ag)/2 


9x(is(A2  -  Ax))  =  e  9  (s)  . 


.  ,  VtT  Z'-siz 


=  1 . 3 . 


. 


1  '  ’  - 


-  Ik  - 


Thus 

is(Ax 

9  (s)  =  e 

r 

(2.18) 

ias  (A 

9p(,as )  =  e 

cpx(is(A2  -  Ax))  ; 


9x(^s(A2  -  Ax)) 


Applying  (2.18)  to  (2.16)  and  dropping  the  subscript  x  on  cp  (0)  , 

X 

Jie(A  +A  )  ie[A. ( 1-a. )+|a- (A  +A  ) ] 

e  9  (ie(A2-Ax))  =  Jt1  e  112  ) ) 


+  «2  6 


ie[A  'l-a0)+|a  (A.+A  )] 

2  1  2  cp  (Ja20(A  -At))  , 


or, 


_0[A  (1-a. )-J(l-a1 )(A  +A  )] 

9(J0(A2-A1jj  =  Jr1  e  cptla^A^) ) 


+  *2  e 


i@[A  (l-a_)-i(l-a  )(A  +A  )] 

2  2  1  1  2  9(ia2e(X2^1)) 


Put  u  =  i  (A2-A1)9 


(2.19)  9(u)  =  it. 


■iu(  l*a^ )  iu(l-a2) 

cp(axu)  +  n0  e  cp(a-2u) 


Because  (2„19)  is  completely  independent  of  and  Ag  ,  i.:  is  obvious  that 

these  parameters  merely  served  to  define  a  scale  factor.  The  change  of 
variable  (2.I7)  transforms  the  old  boundary  conditions  to 


(2.20) 


F(x)  =  0  ,  x  <  -1  ; 

F(x)  =  1  ,  x  >  1 


A  distribution  function  may  be  obtained  by  applying  the  fundamental 
inversion  theorem  to  its  charactei  istic  function  (c.f.  Kendall  [4]  p.9^)« 


|u2-P)9 
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This  inversion  theorem  states  that 


(2.21) 

F(x)  -  F(0)  =  1/(2*)  /  Wt)(l  -  e'lxt)/  (it))  dt  , 

J  ~oo 

where  F(x) 

is  the  required  distribution  function  corresponding  to  the 

characteristic  function  Cp( t ) . 

Inverting  equation  (2.19)  by  the  formula  of  (2.21),  the  left 
member  becomes 


(2.22) 

00 

L  (9(t)(1-e'iXt)/(it))  dt  =  F(x)  -  F(0)  ; 

the  first  term  of  the  right  member  becomes 

it,  r°°  -it(l-a1)  . 

2 ~J  e  (l-e"1X  )/( it ) )  dt 

which  may  be  written  as 


rt,  ,  -00  -it ( 1-a.  )  -it (x+l-a1 ) 

2^  J  {cP(a.1t)[-(l-e  )  +  (l-e  )]/(it))  dt  , 

or,  putting 

u  =  a^t 

(2.23) 

tc,  -iu(l-a1)/a.  -iu(x+l-a1  )/a1 

^  {q>(u)[-(l-e  1  +  (l-e  1  1)]/(iu)}du 

•J  mmOO 

=  ^[-F^l-a^/c^)  +  F(0)  +  F((x+l-a1)/a1)-F(0)] 

=  Jt1[F((x+l-a1)/a1)  -  F^l-a^/a^)]  ; 

while  the  second  term  of  the  right  member  becomes 


•  n  re.  :  5,-  ;  •  sirf'C 


'  "  -  '■  :  t"" 
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*2  o00 

2* 


f  (<p(a2t)  eit^1"a2^(l-e-lxt)/(lt) }  dt 
j  -00 


Jt2  ,°° 

2it 


/  {cp(a  2t)[-(l-« 

J  -00 


-it(ap  -1)  -it(x-l+a0) 

)  +  (l-e  2 


)]/(it))  dt  , 


or,  letting  v  =  a^t 

r00  -iv(ou  -l)/a  -iv(x-l-Kx  )/ap 

(2.24)  gf/  {<p(v)[-(l-e  ^  2)  +  (l-e  2  2)  ]/(lv) )  dv 


-00 


=  Kp[F((x-l4a2)/a2)  -  F(ap-l)/ap)] 


2  '  2 ' 


Collecting  the  inverted  terms  (2.21),  (2.22),  (2.23)  and  (2.24)  the  result 
is 

F(x)-F(0)  =  «1[F((x+l-a1)/a1)-F((l-a1)/a1)]+rt2[F((x-l-Hx2)/a2) 

-F((a2-l)/a2)] 

Putting  x  =  -1  and  imposing  the  boundary  condition  (2.20),  we  obtain 


-F(0)  =  -jr][F((l-a1)//a1)  -  ^((ag-l)/^)  . 


Thus  the  transformed  limiting  distribution  function  satisfies  the  homogeneous 
functional  equation 

(2.25)  F(x)  =  jt1F((x+l-a1)/a1)+jt2F((x-l4a2)/a2)  . 

The  general  solutions  for  the  limiting  distribution  functions 
are  as  yet  unknown,  but  the  properties  demonstrated  by  the  few  published 
results  indicate  that  these  functions  exhibit  rather  unusual  behaviour. 
Because  most  of  these  results  were  derived  directly  from  the  characteristic 


'•)  -  V 

'  .  •. )  ixrsc 
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functional  equation  (2.16)  or  from  the  properties  of  the  operators 
and  ,  explicit  distribution  functions  were  found  in  only  a  very  few 
cases.  Furthermore,  the  known  results  are  generally  restricted  to  equal 
alphas  (a  =  or  equal  pis  (t^  =  =  i)  . 


The  following  review  outlines  the  major  contributions  to  the 
present  knowledge  of  the  properties  of  the  limiting  distributions  of  response 
probabilities.  Jessen  and  Wintner  [5]  proved  that  in  the  case  of  equal 
alphas  and  equal  pis,  the  distribution  functions  are  either  purely  singular 
or  absolutely  continuous.  Kershner  and  Wintner  [6]  showed  that  for  0  <  a  < 
the  limiting  distributions  are  purely  singular  and  their  frequency  functions 
F(x)]  are  non-zero  only  on  a  set  of  points  of  Lebesgue  measure  zero. 

In  particular,  the  frequency  function  for  the  case  a  =  -^  is  non-zero  only 
on  the  Cantor  Ternary  Set.  They  also  proved  that  for  a  =  ^  and  , 

the  distribution  function  is  absolutely  continuous  and  indeed,  in  this  case, 
F(x)  =  i(x+l)  ,  (-1  <  x  <  l)  .  Additional  properties  for  the  case  a  =  ± 
were  derived  by  Kemeny  and  Snell  [7].  These  authors  proved  that  for  all 
^1  ^  n2  ’  fc^e  distributions  are  again  purely  singular.  Karlin  [3] 

generalized  many  of  the  above  results.  He  was  able  to  prove  that  the  limiting 
distributions  are  either  purely  singular  or  absolutely  continuous  for  all 
possible  values  of  the  parameters.  Furthermore  these  distributions  are 
purely  singular  for  all  <  1  . 

For  the  remaining  cases  of  >  1  ,  less  is  known.  Wintner 

[8]  showed  that  for  the  case  of  equal  pis  and  alphas  with  a  =  2 


(n  =  1,2,3,  . ..),  the  limiting  distribution  functions  are  absolutely  continuous 


< 

.  .  woo  ■  ot 

'  -V  | 
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and  even  possess  n-  1  continuous  derivatives.  McGregor  and  Hui  [2] 
demonstrated  that  in  the  equal  pi  and  alpha  case,  for  a  =  1//2  ,  the 
characteristic  functional  equation  (2.19)  could  be  written  as  the  character¬ 
istic  function  of  a  convolution  of  uniform  distributions  over  unequal 
intervals.  Using  this  result,  the  distribution  function  was  found  to  be 
composed  of  three  piecewise  continuous  polynomials  of  degree  two,  Zidek  [9] 
generalized  this  approach  and  showed  that  for  rt^=jtg=|-  ,  and  a  =  2 
(n  =  3>4,5,  •••)  >  the  distribution  functions  consist  of  2n-l  piecewise 
continuous  polynomials  of  degree  n.  Finally,  in  an  early  paper,  Erdos  [10] 
proved  that  for  jc^=jr2=i  >  there  exists  a  set  of  alphas  in  the  interval 
i  <  a  <  1  which  yield  purely  singular  distributions.  These  values  of  a 
lie  between  the  values  which  yield  the  absolutely  continuous  piecewise 
polynomials . 

§2.4  Some  Additional  Properties  of  the  Limiting  Distributions 

In  the  previous  section  it  was  shown  that  the  distribution  functions 
of  the  limiting  distributions  of  the  transformed  response  probabilities 
satisfy  the  functional  system* 

(2.26)  F(x)  =  jc1F((x+l-a1)/a1)-Ht2F((x-l4a2)/a2)  ,“ 

(2.27)  F(x)  =0  ,  x  <  -1  ;  F(x)  =  1  ,  x  >  1  ; 

(2.28)  Jtj+rtg  =1  ,  0  <  rtprtg  <1  ,  0  <  a1  ,  a2  <  1  . 


* 


Functional  equation  and  auxiliary  conditions. 


. 
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Theorem  2.1  For  <  1  the  solutions  of  the  functional 

system  (2.26)  -  (2.28)  are  purely  singular  functions  which  are  constant 
everywhere  except  at  a  non -denumerable  set  of  points  of  Lebesgue  measure 
zero . 

Proofs  The  proof  of  this  theorem  is  composed  of  three  distinct  parts. 

In  the  first  part  an  algorithm  for  constructing  the  exact  distribution 
functions  is  derived.  This  algorithm  shows  that  the  distribution  functions 
take  constant  values  on  a  set  of  closed  intervals  lying  in  the  domain  (-1,  l) 
The  second  part  of  the  proof  comprises  Lemma  2.1  which  shows  that  the  set 
of  intervals  defined  by  the  algorithm  are  disjoint  for  all  <  1  . 

The  third  part  of  the  proof  starts  with  Lemma  2.2.  In  this  lemma  it  is 
proved  that  the  set  of  disjoint  intervals  covers  the  entire  domain  (-1,  1) 
except  for  a  non-denumerable  set  of  points  of  Lebesgue  measure  zero.  Finally 
the  proof  of  the  theorem  (p.  27  )  follows  from  the  fact  that  the  distribution 
can  be  uniquely  constructed. 

In  the  open  interval  -1  <  x  <  1  ,  the  arguments  of  (2.26)  are 
ordered  such  that  (x+l-a^)/a^  >  x  >  (x-l+a^)/^ .  Since 


(x+l-a^/a^  >  1 

when 

x  >  2a, 

and 

(x-l-t<x2)/a2  <-l 

when 

IA 

h-1 

the  boundary  conditions  (2.29)  define  a  solution  if  an  interval  exists 
such  that 

1  -  2cg2  >  2<x^  -  1 


■ 

. 

1  '  ;  r_  ,  .  , ' :  :  '  ■'  • 
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This  condition  is  satisfied  if  and  only  if 

al  +  a2  <  1 

Thus  for  any  x  in  the  interval 


(2.29) 

2a^  -  1  <  x  <  1  -  2a2 

we  have 


(2.30) 

F(x)  =  x1  ,  [2a1  -  1,  1  -  2ag]  . 

From  the  above  solution  function,  two  new  intervals  and  their 
corresponding  solutions  can  be  constructed.  One  of  these  intervals  is 


defined  by 

2a ^  -  1  <  (x+l-a^)/a^  <  1  -  2ag  . 

Solving  for  x 

gives  the  interval 

(2.31) 

a^(2cc^-l)  “  (I"00!)  <  x  <  a-^l-Sa^)  "  (l-a^)  . 

The  other  interval  is  given  by  the  condition 

2a. x  -  1  <  (x-l+a2)/a2  <  1  -  2a2  , 

which  defines  the  interval 

(2.32)  a2(2a1-l)  +  (l-a2)  <  x  <  a2(l-2a2)  +  (l-a2)  . 

The  interval  (2.31)  lies  to  the  left  of  (2.29)  and  is  disjoint 
from  it.  The  contrary  assumption  that 


<x^(l-2<x2)  -  (l-a^)  >  2a^  -  1 
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leads  to  the  contradiction 

-2a^a^  >  0  ,  since  (0  <  ,  a2  <  l) 

Similarly  the  interval  (2.32)  lies  to  the  right  of  (2.29)and  is  disjoint 
since  the  contrary  assumption 

1  -  2a2  >  a^oc^-l)  +  (1  -  ag) 

yields 

0  >  2a^a2 

Thus  the  new  intervals  are  disjoint  for  all  acceptable  values  of  and 

<x„ . 

2 

The  solutions  over  the  new  intervals  are  defined  as  follows; 

i)  when  x  e  (2.31)  ,  (x-l4a2)/<x2  <  -1.  Substituting  into  the 
functional  equation  (2.2 6)  gives 

(2.33)  F(x)  =  .rtj2  ,  [a1(2a1-l)  -  (l-a^),  ^(1-2^)  -  (l-a^)]  . 

ii)  when  x  e  (2.32)  ,  (x+l-a^)/a^  >1  so  the  functional  equation 

gives 

(2.34)  F(x)  =  Jf1  +  Jtgitj  ,  [ag(2a1-l)  +  (l-ag),  ag(l-2a2)  +  (l-a2)] 

In  general,  if  xg  is  an  endpoint  of  an  interval,  new  endpoints 

x  ’  are  determined  from  x  by 
e  e 

(x  ,+l-a1)/a1  =  x  giving  x  '  =  a,x  -(l-a.)  , 

e  1  '  1  e  e  lei 

(xe'-l-wo2)/a2  =  xg  giving  xg '  =  a2xe+(l-a2)  . 
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For  later  convenience,  define  the  two  transformations 

(2.35)  BjX  =  oOjX  -  (l-ax)  , 

(2.36)  B2x  =  agx  +  (l-a2)  . 

At  any  given  stage,  the  operators  B.^  and  B2  generate  the  new  intervals 
by  a  mapping  of  the  intervals  determined  at  the  previous  stage.  In  particular, 
the  new  endpoints  are  given  by  the  transformations  B^xg  and  BgX^  ’  where 
{x^}  represents  known  endpoints. 


Applying  the  algorithm  a  third  time  gives  four  new  intervals  with 
their  corresponding  solutions. 


(2.37) 

F(x)=jt15  , 

[a12(2a1-l)-a1(l-a1)-(l“a1),a12(l-2a2)“a1(l-a1) 

-(l-a1) ]  ; 

(2.38) 

\  2  2 

Ftxj^rt^  it2  , 

[a1a2(2a1-l)+a1(l-a2)-(l-a1),a1a2(l-2a2)-Kr1(l-a2) 

-(l-a2)]  ; 

(2.39) 

2 

F(x)=Jr1+it2Jt1  , 

[<X2a1(2a1"l)  "a2  ( 1“^  )+(l-a2 )  ,0^0^  (1‘2cl2  ^  "a2  ^  1"al ) 

+  (l“a2^  > 

(2.1t-0) 

2 

[a22  (2^-1  )-kx?  ( l-a2 )+( l-a2 )  ,a22  ( l-2a2  )-hx2(  l-a2 ) 

+(l“0t>2)  ]  . 

th 

Since  each  interval  at  the  n  stage  is  defined  by  the  trans¬ 
formation  B^  or  B2  applied  to  the  endpoints  of  an  interval  found  at 
the  n-lSt  stage,  there  must  be  2n  ^  intervals  to  be  found  at  the  nth 
application  of  the  algorithm.  This  procedure  for  constructing  the  solution 
may  be  repeated  as  often  as  required  to  outline  the  distribution  function  to 


■ 

- 

■ 
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the  desired  degree  of  accuracy,  provided  only  that  the  intervals  do  not 


overlap. 

Lemma  2.1: 

The  intervals  generated  in  the  above  algorithm  for 

constructing 

the  solutions  of  the  functional  system  (2.26)-(2„28)  for 

a^-f<x2<l  are 

disjoint . 

Proof  s 

The  following  properties  of  the  transformations 

and  are  required  for  the  proof. 

i)  B^x  =  a^x  ~(l-a^)  and  B^x  =  <x0x+(l-a,-,)  are  linear  functions 
of  x  .  Thus  B^x*  >  B^x  and  B^x'  >  B^x  if  and  only  if  x*  >  x  . 


ii) 

B.,x  <  x  for  all  x  e  (-1,  1)  ,  since  by  assuming  the  contrary 

that 

<x.^x  -  (1-a^)  >  x 

gives  x  <  -1  ,  which  is  a  contradiction.  Furthermore,  B^x  e  (-1,  2a^~l) 


for  all  x  e 

(-i,  i). 

iii) 

B^x  >  x  for  all  x  e  (-1,  l)  .  Again  the  contrary  assumption 

a£x  +  (l-a2)  <  x 

gives  the  contradiction 

x  >  1 

In  this  case  B^x  €  (l-2<Xg,  l)  for  all  x  e  (-1,  1)  . 

Given  an  interval  (x,y)  =  I  ,  define  to  be  the  interval 

(B^x,  B^y)  and  B^I  to  be  the  interval  (B^x,  B2y)  . 


.  '  i  ,  :  '■  ■  -  t'l  ■'  <  ■  ,  3.  {  1  £  1 
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At  the  first  stage  in  constructing  the  solution,  the  functional 
equation  and  boundary  conditions  define  the  interval  (2<x^-l,  l-2a^)  . 
Denote  this  interval  by  I-  -  .  The  remaining  intervals  at  this  stage, 
namely  (-1,  2a^-l)  and  (l-2ttp,  l)  we  shall  call  gaps  and  denote  them 
by  G1(1  and  G^. 

At  the  second  stage,  consider  the  effects  of  the  transformations 
and  .  B^  maps  the  domain  (-1,  l)  onto  (-1,  2a^-l)„  Thus 

G1,1^B1!G1,1^I1,1UG1,21  •  Denote 


B1G1,1 

by 

°2,1 

} 

B1I1, 1 

by 

¥,1 

> 

and 

B1G1,2 

by 

G2,2 

9 

Since  B^  maps 

Denote 

(-1,  l)  onto  (1-2^, 

1)  we 

have  G1;2  Db2[G1>iUI1>1  UG1j2] 

B2G1,1 

by 

G2,3 

f 

¥i,i 

by 

X2,2 

> 

and 

B2G1,2 

by 

G2,4 

• 

Thus  I  .  and 

I  are  bounded  by 

2,2 

G2,l(1= 

1,2, 5, *0  • 

The  same  argument 

also 

holds 

at  the 

third  stage.  Since  B^G^  ^ 

maps  onto  ^  then  B^[Gp 

bi[G2,3vjI2,2UG2,4]  must 

,iV 

map 

G2,2^ 
onto  G2)2  . 

must  map  onto  G  .  .  Also 

^  y  J- 

Thus  the  new  intervals  and 

gaps  generated  by  B^  at  the  third  stage  may  be  denoted  as  follows. 
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B1G2,1 

by 

°5,1 

y 

B1I2,1 

By 

V 

y 

B1G2,2 

by 

V 

y 

B1G2,3 

by 

G3,3 

y 

*1*2,2 

By 

r3.2 

y 

and 

B1G2,* 

by 

°3^ 

• 

The  new  intervals 

generated  by 

B2  are 

located 

such  that 

G2,3  ^ 

B2*-G2,lW 

1  UG  2, 

2i  , 

G2  4  2-  B2[G2  3U  T2  2U  G2  4] 

Continuing  the  notation,  denote 

B2G2,1 

by 

y 

B2^2, 1 

by 

X3,3 

y 

B2G2,3 

by 

G3,6 

y 

B  0G0 

2  2,3 

by 

G3,7 

y 

B2*2,2 

by 

V 

y 

and 

B2G2,4 

by 

G3,8 

o 

From  this  scheme  it  is  obvious  that  in  the  general  case,  any  new 


interval  will  be  bounded  by  two  gaps 


-  26  - 


Figure  2.  Formation  of  the  Intervals  of  Solution  in  the  first 

Three  Stages 

Lemma  2.2  s  In  the  limit  as  n  -»<»  ,  the  set  of  disjoint  intervals 

1^  ^  (m=l,2,3,  nj  i=l,2,  2m”'*‘)  covers  the  domain  (-1,  1)  except 

for  a  non-denumerable  set  of  points  of  Lebesgue  measure  zero« 

Proof «  Given  an  Interval  [xp  x^]  of  length  £  =  Xg-x^  , 

maps  the  interval  onto  the  new  interval  [a^x^"(l-a^),  a^x^-Cl-a^) ] 

length  a1(x2-x1)  =  c^i  .  Similarly  Bg  maps  the  given  interval  onto 

[a2x1+(l-a2),  a2x2+(l-a2)]  of  length  a2(xg-x^)  =  .  Since  both  B^ 

and  B^  map  every  interval  I  ...  to  define  the  I  .  ,  the  sums  of  the 

2  r  J  m-1,1  m, j 

1 

lengths  of  Im  (j=l,2, 2  “  )  is  a^-w-g  times  the  sums  of  the  lengths 

m  o 

of  I  ,  .  (i=l,2, . . .,2  ™  )  .  Thus  the  sums  of  the  lengths  of  [I  .  )  over 
m- 1 ,  i  m, i 
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m  gives  a  geometric  series  with  initial  value  (l~2a  )-(2a1~-l)  =  2(l-a.-a0) 

C.  1  JL  C. 

and  common  ratios  (a^-KX0).  Summing  the  series 
m 

(2.41)  lim  V  2(l-a1-a2)(a1-KX-2)n  =  2(l~a1~a2)/[l-(a;L4a2)]  =  2  . 

m->  oo  l— j 
n=o 

This  value  is  the  sum  of  the  Lebesgue  measures  of  all  the  intervals „ 
Tlius  the  sequence  of  disjoint  intervals  covers  (-1,  l)  except  for  the  set 
of  points  of  Lebesgue  measure  zero  comprising  the  gaps  of  the  limiting  case. 

The  number  of  gaps  at  the  n^  stage  is  2R  .  Thus  in  the  limit, 
the  number  of  gaps  becomes 

lim  2^ 

N— >  00 

In  the  limit  N  is  of  cardinality  yC0  .  A  well  known  theorem  of  trans- 
finite  numbers  states  that 

(2.42)  2X°  =  c 

where  c  is  the  cardinality  of  the  continuum. 

Proof  of  Theorem  2,1 

The  proof  follows  immediately  from  lemmas  2.1  and  2.2  and  by 

noting  that  the  solutions  given  by  the  algorithm  are  constant  over  the 

intervals  fl  . )  . 

n, 

Thus  for  the  case  <x^4a0  <  1  the  limiting  distributions  of  the 
transformed  response  probabilities  are  uniquely  defined  and  can  be  constructed 
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to  any  desired  degree  of  accurracy. 

In  the  case  a, +a_=l  ,  the  intervals  I  ,  of  theorem  2.1  shrink 
12'  n,  i 

to  points,  since  1-20,2  =  2a^-l  .  However  the  algorithm  for  constructing 
the  solution  may  be  used  to  find  the  values  of  the  function  at  a  non- 
denumerable  set  of  points.  In  addition  to  the  construction,  we  have  the 
following  special  result. 

Theorem  2.2:  Necessary  and  sufficient  conditions  that  the  limiting 

distribution  be  the  uniform  distribution  are  ,  cXg  =  ^  t 

"J'jf 2=  -0  • 

Proof  of  Sufficiency. 

The  theorem  will  be  proved  by  showing  that  the  moments  of  the 

distribution  for  a,,=jr.  and  ao_=jt~  are  the  moments  of  the  uniform  distribution. 

1  i  2  2 

Before  the  main  theorem  can  be  proved,  we  again  require  an  intermediate  result. 

Lemma  2.3?  The  characteristic  functional  equation  (2.19) 

-iu(l-a^)  iu(l-ag) 

<P(u)  =  rt-je  cpCa-ju)  +  *2  e  <P(a<2u)  , 

with  9(0)  =  1  defines  a  unique  set  of  moments  Vm  where 

,m 

(Dra  Vm  -  ±-m  [*(«)] 

du 

Proof:  Putting  u  =  0  and  substituting  cp(0)  =  1  In  the  above  functional 

equation  merely  reaffirms  the  identity 

v  —  1  . 


Differentiating  the  functional  equation  with  respect  to  u  gives 


:  v"  '-.s: 
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-iu(l-a  )  -iu( i-a  ) 

cp*(u)  =  a^e  9'  (a][u)  -  a^e  cp'(a2u) 

-iu(l-a  )  iu(l-a  ) 

+  i[-n1(l-0'1)e  cpCa^)  +  it2(l-a2)e  9(a2u)] 

Therefore 

9'(0)  =  i[-jt1(l-a1)  +  jr2(l-a2)  ]/[  l-a^-a^  ]  . 

cp"(o)  may  be  evaluated  by  differentiating  the  functional  equation  twice 

and  substituting  for  cp(o)  and  cp'(O).  Repeating  this  procedure,  all  the 
dm 

terms  — -  cp(u)l  (m=l,2. 3* . • • )  can  be  evaluated  in  turn.  Thus  the 
,  m  v  u=o  1  1  7 

du 

moments  are  uniquely  defined  by  (2.19). 

Proof  of  Sufficiency; 

The  theorem  will  be  proved  by  showing  that  the  moments  of  the 
distribution  for  cx^  =  and  =  xc2  are  the  moments  of  the  uniform 
distribution. 


The  characteristic  function  of  the  uniform  distribution  over 
(-1,  1)  is 

,1 


9(t 


>  -  if 


XXf 

e  dt 


=  i  [eU  -  e“U]/(it) 


(2.43) 


cp(t)  = 


sin 


M 


Thus  the  moments  V  are  given  by 

m 


■  I  ,  '  J  .  : 


30 


(2.44) 


v  =  [1  +  (-l)m]/[2(m+l)]  . 


m 


Bush  and  Mosteller  ([1]  p.88)  give  the  following  recurrence 

formula  for  the  moments  of  the  distribution  of  response  probabilities  at 
s  t 

the  n+1  trial  as 


m 


m 


„  V  /m\  m“u  uw  V  /m\  m“u  UW 

V  i  >  (uja-  a.  V  +jt0  >  (u)a^  a.  V 

m,  n+1  l^x/l  1  u,n  2  ^  2  2  u,n 


u=o 


u=o 


In  the  limiting  case,  the  recurrence  formula  becomes 


m 


m 


V  m_U  UT7  V  /Hlx 

Vm  “  (u)  *1  “l  Vm  +  *2  ^  (u) 


mx  m-u  uTI 

a0  a 0  V  , 

2  2  m  ’ 


u=o 


u=o 


where  V  =  lim  V 

m  m,n 

n— >  oo 


Solving  for  V  gives 
m 


m-l 


m-1 


„  /,  m  nu  ST1  ,m,  m-u  u„  V"*  /m\ 

Vm(1-nlal  -«2a2  >  =  *1^  (U)  31  “l  Vu  +  "2  2,  (U) 


m-u  u 
a2  a2  *u 


u=o 


u=o 


Transforming  the  above  equation  so  that  it  defines  the  moments 
for  the  transformed  distribution  over  (-1,  l)  and  putting 


oc^  =  a 


a2  =  1  -  a  , 


rt]_  -  a  ) 


«2  =  1  -  a-  > 


yields 


m-l 


m-l 


(2.45)  Vm(l-am+1-(l-a)m+1)  =  (u)  (-l)m”UaUVu  +  (l-a)  ^  (u)a,m  U(l-a)UVu 


u=o 


u=o 


o  ---  ■ 


-O  ■  VC -■  ■  1  >.  .V-  ..  ~  '  ~  ,  "  .■  :  :■  ■  =  u v 
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By  definition  V  =  1  . 

o 

The  induction  hypothesis  is  =  [l+(-l)m]/[2(m+l)].  For 
m=l  ,  (2.1+5)  gives 


V1[l-a2-(l-a)2]  =  -a(l-a)V  +  (l-a)aV  =  0  . 


Therefore 


vx  =  0  . 


Assuming  the  induction  hypothesis  to  hold  for  n  <  m-1  ,  (2.1+5)  yields 

m-1 

TT  r-  m+1  /.  \ m+1 i  \  /tn>.  ,  .  \tn-U/-  \m-u  u+lr_  f  „  iUi  „  • 

Vm[l-a  -(1-a)  ]  =  i  )  (u)(-l)  (1-a)  aL  [i+(-i) 


u=o 

n-1 


i  Y  (S)am-u(l-a)u(l+(-DU]/(u+l)  , 


u=o 


m-1 


1  V  /m+l\/  \tn-U  u+l/i  /  n\u\ 

-  2+sry  1  (u+i)(-1}  (1-a)  a  (1+(*1)  ) 


u=o 

m-1 


1  V-1  ,m+l*  m-u,.  \U+] r  ,  .  xu. 

+  2+STy  L  ^  a  (  a)  [1+(_1)  1  * 


u=o 


Thus 

(2.1+6) 


m 


Vm[l-am+1-(l«x)m+1]  J  (l+(-l)v-1)(ml1)[(.l)m+1-v(l-a.) 

v=i 

X 


m+1  -  v 


v  m+1  -v,.  <v, 

a  -kx  (1-a)  J 


Case  I  :  m  =  2k+l  (k=0, 1,2, . . . )  . 


i)  When  v  is  even,  the  corresponding  term  of  the  sum  vanishes. 


ii)  When  v  is  odd,  say  v  =  2£  +1  ,  (.£=0,1,2,...),  the  summand  is 


'  a  ■'  •■!  bqJLI  ■  ; 


'■  m  ■:  v  ■  ,r.-  ■  . .  v 


•  ' ...  :  ■■  r  .  r  !  ,  \ 


I+'IK 


c  -  J 

■  u«;  e  ' 

r  >o 
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,2k+2u  v2(k-i)+l  2UI  2(k-£)+l,.  v2i+l, 

2i+l  ”  (1_a)  a  +a  (1-a)  ] 


Now  let  v  =  2(k+l)-2i-l  =  2(k-i!)+l  .  The  summand  becomes 

/  2k+2  \  r  /,  n2X+1  2(k-i)+l  2.0+1, ,  x2(k-i)+l, 

^(k+D-ei-iJ  [-(1‘a)  a  **  (1‘a)  ] 

In  this  case  the  terms  of  the  sum  cancel  out  in  pairs  so  V„,  ,=0  . 

2k+l 

Case  II  :  m  =  2k  (k=0, 1,2, . . . )  . 


i)  When  v  is  even,  the  corresponding  term  of  the  sum  vanishes. 


ii)  When  v  is  odd,  say  v  =  2.0+1  (i=0, 1,2, . . . )  ,  equation  (2.46) 


becomes 


k-1 


„  ri  2k+l  x2k+l,  1  V  /2 k+lu,.  x2(k-i)  2,0+1  2(k-i) 

V2k[l-a  -(1-a)  \  (£i+1)t(l-a)  'a 

£= o 


/ ,  >.2,0+1, 
x  (1-cc)  J  . 


The  sum  on  the  right  hand  side  is  simply  the  binomial  expansion 


, »  \2k+l  2k+l  /-  \2k+l  -  2k+l  ,,  \2k+l 

( 1-a-ta)  -  a  -  (1-a)  =  1-a  -  (1-a) 


Therefore 


2k  2k+l 


and  sufficiency  is  proved. 
Proof  of  Necessity> 


The  zeroth,  first  and  second  moments  impose  the  conditions 


'  '  '  '  ■  vf 


sq  y  o  sn:  yJ  g rsikr  ■  sc  •  ■  .  rs 


'  1  ;  <>  Q  : 


a/  e  \  I  r  •  ,  : 


, 
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(2.1+7) 

*1  +  *2  =  1 

(2.48) 

^1(1‘a1)  =  *2(l-a2) 

(2.49) 

3[jt1(l-a1)2  +  it2(l-a2)2]  = 

From  (2.48) 

(2.50) 

(l-a2)  =  it^l-c^)/^  , 

whence 

(2.51) 

a2  =  1  ‘  Jtl(1_ai  Vrt2  * 

Substituting  for 

a2  ,  l-a2  and  *2  in  (2.49) 

3[3t1(l-a1)2+(l-Ir1)Jt12(l-a1)2/(l-Jt1)2]  =  l-*^2 

2  2 

3jr1(l-a1)  [l+jr^l-j^)]  =  1-Jt1a1  -  (l-3t1 )  +  2^1(l-a1)  -  * 

=  Jt1(l-a1)[3  +  ax  -  jt1(l-a1)/(l-it1 

If  /  1 

3(l-a1)[l+jt1/(l-jt1) ]  =  3  +  ax  -  1^(1 


(l-®!) if 

12(l-a1)2/( 1-Sl) 
)]  . 

"a^  )/(^-”^^)5 


Therefore, 
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From  (2.51) 


Therefore 


This  completes  the  proof  since  we  have  already  shown  the  condition  <x^  =  , 

a2  =  rt2  suf^^c^ent  to  yield  the  uniform  distribution. 

We  now  consider  the  distributions  for  >  1  .  Under 

this  restriction  on  the  alphas,  Lhe  functional  system  (2.26)  -  (2.28)  no 
longer  lends  itself  to  a  unique  construction  of  the  limiting  distribution 
functions.  However,  the  following  theorem  shows  that  the  functions  satisfying 
the  functional  system  are  unique  distribution  functions. 


Theorem  2.3  The  functional  system  (2.26)  -  (2.28)  uniquely  defines 

the  distribution  function  corresponding  to  the  characteristic  functions  which 
sat-isfy  equation  (2,19). 


Proof 1  Lemma  2.3  shows  that  the  moments  generated  by  the 

characteristic  functional  equation  (2,19)  are  unique  so  the  characteristic 
functions  satisfying  (2.19)  are  also  unique.  Furthermore,  the  inversion 
of  the  characteristic  functional  equation  yields  the  unique  functional 
equation  (2.26). 

Suppose  there  are  two  functions  F(x)  and  G(x)  ,  both  of  which 
satisfy  the  functional  system.  Then  the  moments  corresponding  to  these 
two  functions  must  be  the  moments  defined  by  (2.19).  Since  F(x)  and 


. 


rjiinri  -I :■  ::::  i  i  r-oi  tfsm  r.< 

si  j  aworfa  m  10.  -la  g.„  o.r  fo‘. 

.  /TO;:.? 01  .  o  3ir$kj3$lb  a  p  ■  -  ■  '/■/-.  0 o.l  1  •  •  „•  £1 

..  ..  ........ 


. 


v  b  >  .  i  ■  -  ■ . 


-  35  - 

G(x)  have  identical  momenta,  F(x)  =  G(x)  by  the  uniqueness  theorem  of 
characteristic  functions.  ([Ill  p.35)« 

Theorem  2.4  For  >  1  there  exists  an  absolutely  continuous 

function  F(x)  which  satisfies  the  functional  equation 

F(x)  =  ^((x+l-a^/a^)  +  *2F(  (x-l-Kc2)/a2) 

for  all  x  in  (-00,  2a^-l)  and  also  satisfies  the  boundary  conditions 

F(x)  =0,  x  <  -1  ;  F(x>  =  1  . 

This  function  is  piecewise  in  the  sense  that  the  functional  form  is  different 

over  each  interval  of  a  set  of  intervals  (b^,  b  ^)  =  3, . . .  ,M) 

covering  (-1,  l).  In  the  interval  I  ,  the  function  takes  the  form 

m 

m 

(2.52)  Fjx)  =  c[(x+l)P  +  y  K.(x-bi)P]  , 

i=2 

where  p  =  In  jt^/ln  ,  [K^)  are  constants  and  c  is  a  normalizing  constant. 

Proof:  In  the  open  interval  (-1,  1)  the  arguments  of  the 

functional  equation  (2.26)  are  ordered  such  that 

(x+l-a^/a^  >  x  >  (x-l+a2)/a2  . 

Making  the  change  of  variable  u  =  (x+l-<x^)/a^  transforms  the  functional 
equation  to 


F^jU-l-xx^)  =  jt^F(u)  +  jt2F(  (a1u-fa1+a2-2)/a2) 


.  .  .  :  £.  -  ■: 

■  '  '■  VC  :V  '  ■ 


' 


A 

.  •-  >  =-  •. 


,  •'  oc)  :<  X  *•'  ,.X>  ,  :•••  ) 
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Re-arranging  the  terms  in  order  of  decreasing  argument 


(2.53) 


it1F(u)  =  Ftcyi-l-Ki^)  -  jt2F((a1u-Kt1-KX2-2)/a2) 


As  u  increases  from  -1  ,  there  exists  an  interval  such  that 
(a^u4a^-ta2-2)//a2  <  -1  .  This  interval  is  (-1  <  u  <  (2-2a2~a^ )/a^ )  ,  and 
we  shall  denote  it  by  1^  .  For  u  e  1^  ,  the  lower  boundary  condition 
reduces  (2.53)  to 


(2.54) 


jt1F(u)  =  F(a^u-l+a^ )  ,  (u  e  Ij) 


In  order  to  meet  the  condition  F(-l)  =  0  ,  we  try  a  solution  of 


the  form 


(2.55) 


F1(u)  =  c(u+l)P  , 


(u  €  Ix) 


Substituting  (2.55)  into  (2.54)  ,  the  left  member  gives 


CJt 


and  the  right  member  yields 


c(a^u  + 


Rewriting  the  above  term  in  the  form 


ca1P(u  +  l)P 


it  is  obvious  the  trial  function  will  satisfy  (2.54)  if 


al  = 


. 


I  .  |  i 


vig  rte-fpraK  srb  ,  (4?, .2)  pdMl  I 


-  ■  ‘i  :a'c9*J  avode.  i.  J:/  .•  i .  •/;>> 
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Thus  we  must  have 

(2.56)  p  =  In  itj/ln  . 

For  u  greater  than  the  upper  boundary  (2-2a.p-a^)/a^  of  1^  , 
there  exists  a  second  interval  Ip  such  that  for  u  e  Ip  ,  both  the 
arguments  (a^u-l-Kt^)  and  (a^u4a^-KXp-2)/a,p  lie  in  1^  .  The  upper 
boundary  of  the  interval  Ip  is  determined  by  the  condition 
ttjU-l-KX^  =  (2-2a2-a1)/a1  .  Thus  Ip  =  [(2-2^-^)/^  ,  (2-2a2-a^)/a^]  . 

The  function  Fp(u)  (u  e  Ip)  can  be  constructed  from  (2.53) 
by  substituting  the  known  function  (2.55)  into  the  right  hand  member  of 
the  functional  equation  (2.53)*  Following  this  procedure  we  have 

7tiF2(u)  =  c(a1P(u+l)P  -«2[  ^a1u-a1-HX2-2)/a2+l]p]  , 

=  c(a1p(u+l)p-jt2(a1/a2>pfu  i-(a1-KX2-2)/a1+a2/a1]p)  . 

Since  <x^P  = 

(2.57)  F2(u)  =  c{(u+l)p-3t2/a2P  [u-(2-2a1-al)/aJ]  P)  ,  (u  €  Ip)  . 

Let  b2  =  (2  -  2a2  -  a1)/a1  and  . 

Thus  (2.57)  is  of  the  form 

Fg(u)  =  c{(u+l)P  +  K2(u-b2)P)  , 

where  bp  is  the  lower  bound  of  Ip  (c.f.  above). 


;;  50  I 
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In  the  third  interval  I. 


complicated.  The  upper  boundary  of 


,  the  situation  becomes  slightly  more 
1^  is  determined  by  either 


i) 

ii) 


the  argument 
the  argument 


a.u  -  1  +  a.  crossing  into  I_ 

i  l  3 

(a^u  +  -2)/a^  crossing  into 


f 


whichever  occurs  at  the  smallest  value  of  u  .  Thus  the  upper  boundary 

of  I  is  given  by 

3 

Min[  (2-2a2-a15)/a15  ,  2(0,^^^ -a^2 -a^)/^2]  . 

The  function  F,(u)  (u  e  I,)  is  of  the  form 

3  3 

F3(u)  =  c{(u+l)P+K2(u-b2)P4K3(u-b3)P) 

where  b  =  (2  -  2ag  -  a^2)/^2  . 

This  may  be  verified  by  substituting  into  (2.53)  using  the  fact 
that  the  argument  (a^u-l-nx^)  e  I2  and  the  argument  (a^u-nx^-Hx2-2)/a2  €  1^. 

The  calculation  of  F^(u)  proceeds  routinely  as  above  unless 
the  two  smaller  arguments  both  cross  boundary  points  simultaneously  to  give 
the  same  value  of  the  upper  boundary  b^  of  1^  ;  that  is,  if 

"z  ~z.  2  2  2 

(2-2a2-a1^)/a1J  =  2(a1-w>2-a1  -a2  -a^)/^  . 


If  this  possibility  occurs,  then  the  point  so  determined  may  or  may  not 
be  a  true  boundary  point.  Calculating  gives 


Fj+(u)  =  c((u+l)P  +  ^(u-(b2-lHa1)/a1)P  +  K^(u- (b^+l-a^  )/a>1  )P 
-  it2/a2p[(u-b2)p  +  K2(u-(a2b2+2-a1-a2)/a1)P]}  , 


■ 

:  .  ....  j  1:,  •  v •"  f.  .00  '•  ;  j 


where 
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(u-(b5+l-a1)/a1)P  =  (u-(a2b2+2-a1-a2)/a1)P  . 

If  F^(u)  does  not  reduce  to  F^(u),  that  is  if  ^  ,  then 

the  point  determined  is  a  true  boundary  and  the  two  terms  K^(u-  (b^+l-a^  )/oc.^  )P  - 
(]t2/a2P)K2(u-(a2b2+2-a^-a2)//a^)P  may  be  collected  into  one  term  and  written 
in  the  form  K^(u-b^)p  .  However,  if  =(n^/(x^>)K^  ,  then  the  two  terms 
will  cancel,  in  which  case  we  call  the  point  generated  a  pseudo-boundary 
point,  and  the  true  value  of  b^  is  given  by  the  value  of  u  specified 

when  the  smallest  argument  crosses  b,  . 

3 

Now  we  wish  to  find  the  function  F  (u)  over  the  interval  I  =(b  ,b  ,  ) 

m  m  m  m+i 

Suppose  a,u-l+a,  el.,  and  (a,  u-kx,  -kx,  -2  )/a,_  el.  ,(i  >  j)  ,  while  u  e  I 
11  l-l  '1  1  2  '  2  j-1  —  *  m 

The  upper  boundary  of  Im  is  given  by 


(2.58)  bm+1  =  Min[bi+l-a1)/a1,  (agb  ^-c^-a^/c^  1  . 

If  the  point  determined  by  (2.58)  is  a  pseudo-boundary  point,  the  true 
boundary  point  will  be  given  by 


bm+l  =Min[(bi+1+l-<*i)/a1  '  (a2bj+l+2-ara2)'/al ] ‘ 

Since  m  intervals  have  been  generated  at  this  stage,  the  arguments  (a^u-l-nx^ ) 
and  (a^u-Ki,^+a2-2)/a2  must  together  have  crossed  m+v+2|a  boundaries;  where 
v  is  the  number  of  simultaneous  crossing  of  two  arguments  which  give  true 
boundary  points  and  p  is  the  number  of  pseudo-boundary  points  generated. 

Thus 


m=i+j~2-v-2p. 


. 


d  js&.ayox.  issiliMKs  a.;-*  asrfw 


te 


An  Induction  argument  may  now  be  used  to  prove  that  F  (n)(u  e  X  ) 

m.  '  m ' 


has  the  form 


m 


In  particular 


)  Kk'U"bk)  ^  ' 

(b  <  u  <  b  ) 

'  m  —  m+1 

already  been  proved  for 

m  =  1  and  2  „ 

true  for  all  the  functions  F‘k(u)  (k=I,2j,o 

j-1 

k=2 

(u  e  Ij^1)  l 

i  ■=  1 

+  )  Kk(u-bk)P), 
k=2 

(u  €  Ii=1)  . 

Substituting  into  (2,53)  gives 


JL  w  X 

(2,59)  Fm(u)  =  c{(u+l)P  +  ^(u-Cb^l^)/^)1 

k^2 


rt2/a2Pi(v"b2)P  +  )  Kk(w“(a2bk42°a1-a-2)/a1)P]}  ,  (u  e  Iffi) 


k=£ 


The  (j“2)  values  f  (b^-t-l-a^ )/a^ }  are  the  endpoints  generated  when 

a^u-l-fo^  crosses  into  the  intervals  (l^}  (cf.  (2,58)),  Similarly,,  the 
terms  ({a^b^-h2<^<x^a^)/a^)  are  the  endpoints  which  are  generated  when 
(a^u4a^-Kx,2»2)/a2  crosses  into  the  intervals  fl^)  ,  There  are  only 
i+j-2~v~2p  independent  values,  since  2v  terms  will  collapse  into  v 
terms  and  2p  terms  will  vanish.  Thus  the  right  hand  side  of  (2,59)  can 


be  written  in  the  form 


.  41  - 


F 

m 


(u) 


m 


((u+l)1 


VUV 


k=2 


(u  €  I  ) 
'  nr 


i 


and  the  induction  is  complete. 

To  finish  specifying  the  function  F(x)  ,  two  more  details  are 

required: 

i)  Let  Fq(x)  =  0  ,  x  e  (-«,  -l)  , 

ii)  In  the  final  interval  I  =  (b  ,  l)  ,  evaluate  c  by  the 

condition  F  (l)  =  1  . 
nr 

This  completes  theorem  2.4. 

The  intervals  generated  by  the  above  procedure  may  be  illustrated 


geometrically  as  follows. 


Figure  3*  Geometric  Interpretation  of  the  Method  of  Generating 

the  intervals  I.  , 

1 


The  upper  boundary  of  the  m  interval  is  given  by  the  value 
of  the  argument  (x+l-ot,^  when  either  x  or  (x-l+a^/a^  crosses  a 

a  boundary  point  already  determined.  For  example,  the  upper  boundary  of 
1^  is  given  by  the  value  of  (x+l-a^)/a^  when  the  argument  (x-l+a^)/^ 

intersects  the  boundary  b^=  -1  „  Similarly,  the  upper  boundary  of  1^ 

is  defined  when  x  intersects  b^  .  The  procedure  for  finding  the  boundary 

points  continues  along  in  the  fashion  illustrated  until  the  upper  boundary 
of  I  is  greater  than  or  equal  to  +1  . 


"is  .  ''if.  ■>: 
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Theorem  2.3  For  a^+a^  >  1  there  exists  an  absolutely  continuous 

function  G(x)  which  satisfies  the  equation 

G(x)  =  it1G((x+l-a1)/a1)  +  jc2  G((x-l4a2)/a2) 

for  all  x  in  (l-2a2,  °°)  and  also  satisfies  the  conditions 

G(x)  =  1  ,  x  >  1  ;  G ( —  1 )  =  0  . 

This  function  is  piecewise  in  the  sense  that  the  functional  form  is  different 

over  each  interval  of  a  set  of  intervals  (^n+;[>  d^)  =  (n=l,2, 3# « . .  N) 

covering  (-1,  l).  In  the  interval  the  function  takes  the  form 

n 

(2.60)  Gn(x)  =  1  -  B[ (l-x)q  +  ^  Ki(di~x)q]  , 

i=2 

where  q  =  In  jtg/ln  a2  ,  (K^)  are  constants,  and  B  is  a  normalizing 

constant . 

Proof?  Make  the  change  of  variable  v  =  (x-I-kx^)/^  and 

substitute  into  (2.2 6)  to  get 

G(a2v+l-a2)  =  jt1G((a2v-2-a][-a2)/a1)  +  jtgF(v)  . 

Rewriting  in  order  of  increasing  arguments 

(2.61)  tc2G(v)  =  G(a2v+i-cc2)  -  ^  G( (a2v+2-a][-a2)/a1 )  . 

In  the  interval  [-1,  1]  there  exists  an  interval  such  that 
(a2v+2-a1-a2)/a1  >  1  .  Denote  this  interval  [  (2a1-HX2~2)/a2,  1]  by 
.  Imposing  the  condition  G(v)  =  1  ,  v  >  1  reduces  (2.61)  to 


I 


' 


. 


(2.62) 


y 


(v  €  Jt)  . 
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*2G(v)  =  G(a2v+l-a2)  -  Jt]L 

To  meet  the  condition  G(l)  =  1  ,  try  a  solution  of  the  form 

(2.63)  G1(v)  =  1  -  B(l-v)q  ,  (v  €  Jj)  . 

Substituting  the  above  function  into  (2.62),  the  left  member  is 

*2  "  Jt2B^1"V^q  ’ 
and  the  right  member  becomes 

1  -  B(a2-a2v)4  -  7t1  . 

The  above  term  can  be  written 

JC2  -  Ba2q(l-v)q  . 

Thus  the  trial  function  will  satisfy  (2.62)  if 

CCg  =  Jtg  . 

This  imposes  the  condition 

(2.64)  q  =  In  In  a2  . 

For  v  less  than  the  lower  bound  of  =  (2<x^-ia^-2)/a^  ,  there 
exists  a  second  interval  J2  such  that  for  v  e  J2  ,  both  the  arguments 
(<x2v+l-a2)  and  (<x^v+2~(x^-a^)/a^  are  in  .  The  lox^er  boundary  of  J0 

is  determined  by  the  condition 

a2v+l-a2  =  (2a14<x2-2)/a2  . 

2  2 

Thus  J2  is  the  interval  [  (2a^-}<x^  -2)/a 2  ,  2a^+a2~2]  .  The  function 


;  )0  <5  i 3 limoa  9^1  3$m  %  T  i 


\ 
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defined  on  the  interval  may  be  found  by  substituting  into 

the  right  hand  side  of  (2.6l)  to  get 

*2g2(v)  =  l-Ba2q(l-v)q-jt1[l-B(a2/a1)q((2a1-4<x2-2)/a2-v)q] 


Therefore 

G2(v)  =  1-B[  (l-v)q-(jt1/aiq)((2a1-Ki,2-2)/a2-v)q]  . 

Denote  the  upper  boundary  of  J2  by  d2  and  let  =  -jt^/a^q  .  With 
these  substitutions,  the  functional  form  of  ^(v)  becomes 

G2(v)  =  1-B[ (l-v)q+K2(d2-v)q]  ,  (v  e  Jg)  . 

In  the  general  case,  suppose  that  the  argument  (a2v+l-a2)  e  JL  ^ 
and  the  argument  ( (a2v+2-a^-a2)/a^ )  €  ^  (i  >  j)  while  v  €  .  The 

lower  bound  of  is  given  by 

(2.65)  dn  =  Min[  (di-l+a2)/a2  ,  (ctjd  . 

However,  if  the  value  of  d^  determined  by  (2.65)  turns  out  to  be  a  pseudo- 
boundary  point,  then  the  true  boundary  point  d^  will  be 


dn  =  Min[(di+1-l*x2)/a2,  (a^  ^0.^2)/^]  . 

The  same  induction  argument  used  in  theorem  2.4  may  be  applied 
to  prove  the  form  of  G^  is 


G 

n 


(v) 


1“B[ ( 1- v )q  + 


n 


Ki(di  -  v)q]  ,  (v  e  Jn)  . 


■it  j  i.3  (.to. ^  i<.  &bt-  bttsit  s>rf:. 


j.O  ...  -  ..  .  i 
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G(x)  is  completely  specified  by  the  additional  conditions 

Gq(x)  =  1  ,  (1  <  x  <  oo ) 

and  by  evaluating  B  by  the  boundary  value 

y-D  =  o , 

where  is  the  interval  which  contains  the  point  -1. 

Theorem  2,6  If  the  functions  F(x)  and  G(x)  of  the  previous  two 

theorems  are  identical  over  (-1,  1)  ,  then  we  have  constructed  the  unique 
absolutely  continuous  distribution  function  which  satisfies  the  functional 
system  (2.26)-(2.28) . 

Proof:  That  the  functional  equation  is  satisfied  everywhere 

follows  from  the  construction.  The  proof  then  follows  immediately  from 
theorem  2.3  which  shows  that  the  functional  system  (2.26)  -  (2.28)  defines 
the  unique  distribution  function  corresponding  to  the  unique  characteristic 
function  which  satisfies  (2.19). 

If  the  construction  of  theorems  2.4  and  2.5  fails  for  a  given  set 
of  alphas  and  pis,  then  the  distribution  function  may  have  some  other 
absolutely  continuous  form  than  the  one  constructed,  or  it  may  be  a  purely 
singular  function. 
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CHAPTER  III 

THE  SUBJECT -CONTROLLED  MODEL 


§  Introduction 

This  chapter  contains  a  treatment  of  the  Subject-Controlled 
Model  which  follows  much  the  same  outline  as  Chapter  II.  First  the  basic 
model  is  derived,  then  the  equation  satisfied  by  the  characteristic  functions 
of  the  limiting  distributions  is  found.  We  again  standardize  the  character¬ 
istic  equation  by  mapping  the  limiting  distribution  on  the  domain  (-1,1). 
However,  in  this  case  the  standardizing  procedure  gives  two  transformed 
characteristic  equations  depending  upon  the  relative  sizes  of  and 

Applying  the  fundamental  inversion  theorem  we  find  that  the  transformed 
limiting  distributions  satisfy  one  of  two  functional  integral  equations. 

A  study  of  these  functional  integral  equations  reveals  a  few  of  the 
properties  of  the  limiting  distributions. 

§  3.2  Derivation  of  the  Model 

The  Subject-Controlled  Events  Model  will  be  derived  in  a  similar 
manner  to  the  Experimenter  Model  of  Chapter  II.  For  this  model, also,  we 
shall  consider  a  simple  learning  experiment  in  which  each  trial  consists 
of  presenting  a  subject  with  two  mutually  exclusive  alternatives  A^  and 
A^.  As  before,  the  choice  of  A^  is  called  response  R^  and  the  choice 
of  A^  is  called  R^ .  The  major  difference  is  that  outcome  0^  is 
always  associated  with  A^  and  outcome  0 2  is  always  associated  with  A^o 
Under  this  scheme  the  mutually  exclusive  events  are 
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E,  :  R.  followed  by  0.  , 

(3.1)  11  1 

:  Rg  followed  by  0^  . 


In  this  type  of  experiment  it  is  obvious  that  the  outcome  of  any  trial  is 
wholly  determined  by  the  response  of  the  subject. 


Following  the  notation  introduced  in  Chapter  II,  let  p^  be 

the  probability  of  response  and  =  1  -  pR  be  the  probability  of 

til 

response  R^  at  the  n  trial.  The  response  probabilities  p^  and 

are  transformed  by  the  same  operators  and  that  were  used  in 

the  Experimenter  case,  that  is 


(3.2) 

if  E1 

(3.3) 


Pn+1 

occurred  at  trial 

Pn+1 


Q,  p  =  a.  +  a.p 
In  1  In 


n  :  or 


QoP  =  ao  +  aoP 

2  n  2  2  n 


9 


9 


if  E^  occurred  at  trial  n.  The  parameters  ai*a2,al,<X2  satisfy  the  same 
conditions  as  before.  (c.f.  2.7) 

th 

Since  the  subject  determines  which  event  occurs  at  the  n  stage, 

the  probability  of  applying  the  operators  and  is  no  longer  a  constant 

tc i  or  jtg  .  Now  it  becomes  a  function  of  p^  such  that  the  conditional 

s  t 

probability  of  E^  occurring  on  the  n+1  trial  given  that  E^  occurred 

on  the  n*"*1  trial  is  Q^Pn  >  and  c^e  conditional  probability  of  E^ 

s  t  th 

occurring  on  the  n+1  trial  given  E^  on  the  n  is  QgP^. 


As  in  the  Experimenter  Case,  if  the  subject  has  an  initial  response 
probability  pQ,  the  possible  values  of  p^  throughout  the  first  few  trials 


i’ ;  3  -■  &  ,  - 1 
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are  generated  as  shown  in  figure  4. 


TRIAL_ 0_ 1_ 2_ 3 


Figure  4:  Possible  values  of  the  response  probability  p^  and  the  probability 
of  the  occurrence  of  that  particular  value  of  pn  throughout  the 
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first  three  trials.  The  upper  entry  is  the  value  of  p^  in 
terms  of  the  operators  and  and  the  initial  response 

probability  p  .  The  lower  entry  is  Pr ( p  )  also  in  terms  of 
^l,(^2'Po,CIo'  (ReProduced  from  [1]  p.  79, 

As  in  the  previous  chapter, let  the  2n  possible  values  of  the 

T1 

response  probability  at  the  n  trial  be  denoted  p  (m=l ,2, 3, . . . ,2  , 

m,nv  ’  ’  ’  ’ 

n=l,2,3, . . . ) •  Also  let  the  corresponding  probabilities  of  the  occurrence 

of  p  be 
m,n 


=  pr(p„  J 


Thus  if  Pr(R^)  =  Pm  n  on  trial  n  ,  the  possible  response  probabilities 
for  trial  n+1  are  given  by 


(3.4) 


Model . 


This 


^lPm, n  al+alPm,  n 

Vm,n  =  a2*Wn 

completes  the  derivation 


with  probability  p  P  , 

J  m,  n  m,n7 

with  probability  (l-p  )P 

J  x  rm,ny  m,n 

of  the  basic  Subject-Controlled 


§  3.3  Properties  of  the  Model 

The  varying  probabilities  of  the  application  of  the  operators 
and  have  made  this  model  extremely  difficult  to  analyze.  As  a 

result,  very  few  properties  of  this  model  have  been  determined.  However 
Bush  and  Mosteller  [12]  found  an  equation  satisfied  by  the  moment  generating 
function  of  the  distribution  of  response  probabilities  by  an  approach 


.  r  c  .  .  i  a  ■  'o“  ••  1 V'  '•  -r  ’  'j.‘ 
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similar  to  that  of  the  Experimenter  Case. 


Let  mn(3)  be  the  moment  generating  function  of  the  distribution 


of  p  ^  .  By  definition 

«n+l 


m,  n 


m 


n+1 


(9)  ■  I 

m=l 


e  ra,n+lp 


m,n+l  * 


,n 


V  r  (al'K1lP™,n)e  p 

)  I  3  ’  P  P  +e 

Z_i  L  ra>n  m>n 

m=l 


(a  +a  p  )e  "I 

2  2  m’n  (i-p  )p  1 

m  ”  " 


;p  | 
i,n  m,  nj 


,n 


=  e 


a^e 


r  a,p  e  i 

a_@  ^ — 1 

r  a,p  0 

\  1  m,n  _ 

2  \ 

1  m.n  „ 

e  p  P 

+e  > 

e  9  P 

L  m,  n  m,n 

L 

m,  n 

L  7 

m=x 


m=l 


Bse  y 


p_  p 


m=l 


m,n  m,nj 


Therefore 

(3.5) 


mn+l(e) 


(-)e  1  d 

al 


aR 


a26  d 


5?[mn(a10)]+e  2  -n(a2@)-£-)e  ‘  ^[n^e)] 


Thus  in  the  Subject  Case,  the  moment  generating  function  of  the  distribution 

s  t 

of  response  probabilities  at  the  n+1  trial  is  dependent  not  only  on  all 
the  moments  of  all  the  previous  trials  but  also  upon  the  first  derivatives 
of  all  the  previous  moment  generating  functions. 


Replacing  by  7\^(l-cu)  (i=l,2),  equation  (3*5)  becomes 


iC  fo  ■■  !  •'  '  \ 
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1  ^  “®"l  J 

(3.6)  mn+i(e)  =  (-)e  Jg  [mn(<V)]« 


^2(l.-a2)0 


mn  (otp  9 ) 


1  A2(l-a2)0  d 

(~>  6  ¥[mn(a2®)]- 


a, 


Converting  the  moment -generating  functional  equation  (j.6)  to  a  characteristic 
functional  equation  gives 


(3-7) 


W9’ 


i>l  (1-ot! 

(-i/cc^e 


)©d  iA2(l-a2  0 

de  [<PnSe)]+e  '  '  \^2e) 


+ 


(i/a2)e 


iA2(l-a2)0^ 

d0 


[cpn(o20)] 


y 


where  cp  (0)  =  m  (i@)  . 
n'  '  n'  ' 


Harris  [1  p.99l  has  shown  that  a  limiting  distribution  F(p) 
exists  as  n  — >°o  except  in  the  cases  xdiere  either  or  both  of  A2  =  0 
and  A^  =  1  hold.  In  these  exceptional  cases,  all  of  the  probability 
density  will  be  concentrated  at  p  =  0  or  p  =  1  . 

To  begin  the  investigation  of  the  limiting  distributions,  let 


lim  9(0)  =  9(0) 

n— >  00 


The  characteristic  function  of  the  limiting  distribution  of  response 
probabilities  satisfies  the  functional  differential  equation 


iA  (1-a,  )d,  iA  (l-ap)0 

(3.8)  9(0)  =  (-i/a1)e  ^  [9(^0)  ]+e  9(a20)+(i/a2) 

iA  (l-ap)0, 

X  e  ^  [9(a20)]  . 
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As  Karlin's  trapping  theorem  also  applies  to  this  problem,  the 
boundary  conditions  on  F(p)  are 


(3.9) 

F(p)  =0  ,  p  <  Min(*lf*2)  , 

F(p)  =1  ,  P  >  Max(A1,A2)  . 

Following  the  procedure  of  the  Experimenter  Model  we  will  make  a 

transformation  to  map  the  limiting  distribution  onto  the  interval  [-1,1]. 
However  (3-8)  is  not  symmetric  in  a^,a^  and  ApAg  so  t^e  f°ll°winS 
transformations  must  be  applied. 


(3.10) 

X  =  (^i+>2"2p^^l"A2^  *  for  ^2  >  A1  * 

(3.11) 

x  =  (A2+A1-2p)/(A2-A1)  ,  for  >  Ag,  (7^  /  1,^2  /  0) 

Equation  (2.18)  gives 


iio@(A  +A  ) 

cpp(cc0)  =  e  cPxt|ae(7v2='A1)]  ,  (Ag  >  ^)  • 

Thus  for 

\  >  \ 

j  ,  — id-S(^p+A.  ) 

^[<ppM)]  “  Jo  (•  q>x[^e(VhH)  - 

or 

(3.12) 

iicc0(A  +A  ) 

fg  [9p(cc0)]  =  ilct<VVe  q>xC4afl(^2.\)] 

+  •  le  ■ 

The  characteristic  equation  corresponding  to  transformation  (3.10)  may  be 
derived  by  substituting  (2.18)  and  (3.12)  into  (3*8)  to  give 


' 


'  "  '  '  *  ' 


-  _  -  r  •  '■  :'ro  OS'M../! 
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!#(Ap+?0  ie[A1  (l-a- )+ia.  (Tv  +A 

i0[A  (l-a  )-f-ia  (A  +A  )] 

-  (iAX)e  h 


+  e 


i0[A  (1-a  )+ia  (A  +A  ) ] 

-  |(A2+A1)e  2-221  cpJla^A^)  ] 

i-0  [  -4p(  1 "ap )  +  ■gOtp  ) -1  j 

+  (i/a2>e  "  %  WiV^W1) 


Thus 


ii0(l-a  )(A  -A  ) 

“’x'i^W1  =  £<VVe  9x[4al9(A2-X. 

ii@(l-a  )(A  -A  ) 

-  (Vty*  3e 

ii0(l-a  )(A  -A  ) 

+  [l-HA^)]  e  ^  ^  <Px[^2e(\-\)] 


+  (i/a  )e 


iied^KAg-Aj)  d 


&  WVi”)  • 


Writing  u  =  ^(A^-A^)  and  omitting  the  subscript  x  on  cp  we 


(3.13)  <P(u)  =  i(A2+A1)e 


-iu(l-a^) 


cp(axu)  +  [l-i(A2+A1)  ]e 


iu  ( 1  -ar 


-iu(l-a  ) 

-  [Ji(^2-^1)/a1]e  [cp(axu)] 


+  [|i(*2"\  Va2^e 


iu(l-a2)  d 


dJt<P(a2u)]>  ()'S>\) 


)] 

cpx[Jai@(A2-A1)] 


)] 


have 


cp(a2u) 
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The  characteristic  equation  corresponding  to  transformation 

(3.11)  is  similarly  derived  by  interchanging  A^  and  Ag  in  (2.18)  and 

(3.12)  then  substituting  into  (3-8)  to  get 


M0(A  +A  ) 


2 \”i  (l-a-i  )+§ch  (^1  +^0)  ^ 

e  1  2  q>[ie(\-V]  =  HVVe  q>[4a10(ArA2)] 


i0[A  (1-a  )+  |a  (A  +A  ) ] 

-(i/a^e  1  %  (q>[4a10(A1-A2)]) 

i0[A„(l-a„)+|ap(A  +A  ) j 

-i(A1+A2)e  2  2  2  1  2  cp [i^©  (A^)  ] 

i0[A  (1-a  )+|<x  (A  +A  )] 

+dA2)e  2  2  2  1  2  |g  ((Ptlo^e^-^)]} 


Thus 


-1-10  (l-a- ) (A. -A  ) 

9[j0(ArA2)]  =  i(A1+A2)e  ^  9^0 (A^Ag)  ] 

|i0(l-a1 )(A  -A  ) 

-(i/ax)e  %  {cp[ia10(ArA2)]) 


+[1  -  f(VV]e 


|i0(l-a2)(A2-A1) 


9[|a20 (ArA2)  ] 


|i0(l-a2)(A2-Ai)  ^ 

d0  rrLa-g-x-i  "g- 


+(i/a  )e  ^  M^a-pS^-A^)])  • 


Put  v  =  ^0(A1-A2)  to  yield 

iv(l-a1 ^ 

cp(v)  =  i(A1+A2)e  J  cp(a1v)  +  [  l“J(\+^2)  le 

iv(i-a1)d 


-iv(l-a2) 


(a-gv) 


-  [Ji(A1-A2)/a1]a 
+  [^i(A1-A2)/ag]e 


-iv(l-a2)  d 


^  [<P(aiv)] 

^  [cp(a2v)],  (Ax  >  A2,  A1  ji  1,  A2  jl  0) 


(3.14) 


■  .•] 


. 
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Equations  (5.13)  and  (3.14)  are  satisfied  by  the  characteristic 
functions  of  the  transformed  limiting  distribution  functions.  The  transforma¬ 
tions  (3.10)  and  (3.H)  imply  the  new  boundary  conditions  on  the  distribution 
function  F(x)  to  be 


(3.15) 

and 


F(x)=0,  x  <  -1  ; 

F(x)  =  1  ,  x  >  1  . 


Following  the  procedure  of  Chapter  II,  the  fundamental  inversion 
theorem  (2.22)  is  applied  to  equations  (3*13)  and  (3*14).  Carrying  out 
the  inversion  on  (3. 13)  one  term  at  a  time,  the  left  member  becomes 

.  00 


1_ 

2jt 


J  - 00 


{cp(t)[i-e’1Xt:]/(it)}dt  =  F(x)  -  F(0)  , 


while  the  terms  of  the  right  member  become: 


:  1  [ 1-e  1Xt]/(it))  dt 


)  KVV  h  f  MV*' 

J  —CO 

poo  -it(l-a1)  -it(x+l-a1) 

=  i(^!+^2)  '  fcp(a^t ) [ - (l-e  )+(l-e  l/(it)) 


dt 


Putting  0  =  a^t  gives 


4(W  h  / 


-i0'l-a.  )Ai  -j  @(x+l-a.  )/a. 

J  )+(l-e  )l/(ie)j  d 0 


=  i(^1+^2){-F[(l-a1)/a1]+F(0)+F[(x+l-a1)/a1]  -  F(o) )  , 

=  i(\+*2)  fF^x+1"aiVai^  -  F[(i-a1)6hL])  • 


i  - 


i  ■  ,  j  •  ( . 
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ii)  [l-i(A1+A2)]  —  /  (cp(a2t)  e 


it('J"a'  2)  _ixt- 

~  [ 1-e  ]/(it)}  dt 


•,  r00  it  (l-ap)  -it  (x-lm0) 

=  t1-i(A1+A2)]  ^  /  (cp(a2t)[-(l-e  )+(l-e  2)]/(it))dt, 

-7  -00 


=  [i-i(\+^2)]  {F[(x-l-Ki2)/a2]  -  F[  (a2-l)/a2J )  . 

i)  i  r00  j  -it(l-a  ) 

iii)  -  5?  /  tsir  (<f>St)]e  1  [l-e'1Xt]/(it)}dc  , 


which  after  integrating  by  parts  becomes 

i(A  -A  )  r  -it(l-a, )  -it(x+l-a,  ) 

- 2~  "  j  LcP(a-1t)(e  -e  )/(it>] 


t=  -  00 


itfl-a^)  -it (x+l-a^ ) 


Cp(a1t )  {i[  (l-a^ )e  1  -(x+l-a^e 

.00 

-it(l-a^)  -it(x+l-a^)  2 


]/(it:) 


+  [e 


]/(it  )}  dt 


Since  |cp ( t )  |<  1  ,  the  integrated  t^m  vanishes,  leaving 


1 

2a^  2jt 


(1-c^)  /  [cp (ax t )  [ e 

J  -00 

-it (x+l-a^ ) 


■  it ( 1 -a. )  - i t ( x+1 -a. ) 

-e  ]/( it ) )  dt 


[9(a1t)[  (xe 


-it(l-a  )  -it(x+l-a..) 

)/(it)  +(e  -e  )/t  }  dt ^ 


Ag-A^  ri-a,  n°° 
2a 


a.  rr-a1  j.-a.  / 

—  |-g^-  J  (9(a1t)[-(l-e  )  +  (l 


it(l-a^)  -it(x+l-a^) 


)]/(it) j  dt 


-  h(x)f  , 


-  '  -I  -*■  - 
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V\  I  , 


i|(l-a1){F[(x+l-a1)/a1]  -  F[  (l-a^/c^  ]  ]  -  h(x)j  , 

,  f 00 (  -it (x+l-a,  )  -it(l-a1)  -it (x+1-a. ) 

where  h(x)  =  ~  /  i<p(a1t)[(xe  )/(it)+(e  -e  ) 


/t  ]rdt  . 


Since  h(x)  is  uniformly  convergent  in  x  ,  differentiation 
with  respect  to  x  may  be  carried  out  under  the  integral  sign. 


Thus 


h'W  -s 


q>(o.1t) 


•it (x+l-a^ )  _lt(x+1.a  ) 


it 


xe 


-it (x+l-a^ ) 


le 


dt  , 


-  -it (x+l-a.) 

—  /  cp(ait)(xe  )  dt 


x 

2rt 


r00 

/  cp^t)  e 
J  -00 


-it (x+1 -a^ ) 


dt  . 


Put  6  =  a^t  so  that 


oo  -i@(x+l-ct  )/a 

h'(x)  =  -(|-)  J  "  q>(©)e  d0  . 


For  those  values  of  x  for  which  F(x)  has  a  derivative  f(x)  ,  the  fundamental 
inversion  theorem  gives 


(5.16) 


fM  =  57 


/  x  -ixt  , 
cp(t)e  dt  . 
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Applying  (3. 1.6)  to  h'(x)  gives 


h'(x)  =  -(x/a1)f[(x+l-a1)/a1l  . 


Integrating  to  evaluate  h(x)  we  obtain 


hW  «  -f- 


“  £  tf[  (t+l-a1  )/ol]L ]  dt  +  Kx  , 


where  is  a  constant  of  integration, 


Making  a  change  of  variable 


h(x)  = 


^(x+l-a1)/a1 


4.00 


f  (u)  [a^u-l+a^)]  du  + 


=  (l-a^ )F[ (x+l-a^)/a^ ]  - 


(x+1 

”al  )/al 


uf(u.)du  +  , 


or 


h(x)  =  (l-a1)  F[(x+l-a1)/a1]  -  {( (x+l-a^/a^Ft  (x+l-a^ )/a1  ] } 


(x+l-a^ )/a^ 


F(u)du)  +  Kx  +  IC2  , 


where  Kg  is  another  constant  of  integration, 


The  total  contribution  from  iii)  is 


V\ 


2a, 


(l-a1){F[(x+l-a1)/a1  ]  -  F[  ( 1  -oo1  )/ot-1  ] )  -  (1-a^  )F[  (x+l-a1  )/cl1  ] 


+  (x+l-a^)  F[  (x+l-a^  )/a.j  1-a^ 


(x+l-a^  )/a]L 


F(u)du  +  +  Kg 


>  .5  ~ >  rgotnl 


' 
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V\ 

(x+l-a1)F[  (x+l-a1)/a1]  -  ( 1 )F[  (l-a^/o^] 


a. 


( x+1 -a^ ) /a^ 


F(u)  du  +  +I<2 


iv 


2a~  2jt 


(^T  [9(a2t)]  e 


i t ( 1 -a0 ) 


2/(l-e_ixt)/(it)}dt 


Integrating  by  parts  as  before  yields 


{> 


2a„  2«  l[<P(a2t)(e 


)/(it)] 


t=  -  00 


n00  it(l-a  )  -it (x-l-KLp) 

-  ^(agt)  (i[(i-a2)e  +(x-l-ta2)e  ]/( it ) 

J  -00 


it(l-a  )  -it(x-l-«x  )  0 

[e  -  e  ]/( it  ) )dt 


The  integrated  part  vanishes  leaving 


V\  r1_ao  P00  it(l-a  )  -it (x-l4a0) 

-IqH-sr  Lw^t)[e  -e  2i/(it))dt 


P00  -it(x-l+a0)  it  ( l-aa)  -it  (x-1-kx,q)  0 

+  ~  I  9(a2t)[(xe  2  )/(it)+(e  2  -e  )/t2]dt 


Vh 

2a^ 


■|(l-a2)  {F[  (x-l-wo2)/a2]  -  F[  (a^-l  )/ag ] }  +  g(x)j"  > 


where 


■t  r00  -it(x-l-KX  )  it(l-a  )  -it(x-l-ta  ) 

g(x)  =—  I  cp(a2t)[(xe  )/(it)  +  (e  -e  )/t 


]dt  • 


« 


.«»  )  •'  r„A 


-  61  - 


8'  (x) 


1_ 

2jt 


9(a2t) 


-it(x-l+a0) 
e _ 

it 


- i t ( x - 1 +a^ ) 
xe  - 


-it(x-l  -KV,, 
e 

It 


dt 


y 


-x 

27 


°o  -it  (x-l+ctg) 

/  <p(a-2t)e 

4-00 


dt 


y 


(x/ap) f [ (x-l+a2)/a2 ]  (at  the  points  where  f  exists) 


Therefore 


g(x)  = 


g(x) 


i  rx 

—  /  uf [ (u-l+a0)/a0]  du  +  K  , 
a  loo  d 

r  (x-l-Kx2)/a2 


(a2v+l-a2)f(v)  dv  +  , 


=  Kz  -  (l-a0)F[(x-l-Kx2)/a2]  -  ag 


(x-l+a  )/i 


2  2 


3 


vf(v)  dv  , 


-  (l-a2)  F[(x-l4a2)/a2]  -  a2^((x-l+a2)/a2)F[  (x-l-Hx2)/a2] 

^(x-l4a2)/a2 


F(v)  dvr  + 


The  inverted  form  of  iv)  is 


V2 

2a„ 


{(l-a2)[F[(x-l-Rt2)/a2]-F[(a2-l)/a2]]+  K^- ( l-a2  )F[  (x-l-togVogj 

n  (x-l-kx2)/a2 

-(x-l*x2)F[(x-l4a2)/a2]  +  a2  j  F(v)  dx)  , 


|~^  {-(x-l-Kx2)F[(x-l4a2)/d2l  -  (l-a2)F[(a2-l)/a2] 


+ 


F(v)dv  + 


+ 


(where  K-,  and  are 

constants  of  integration). 
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Collecting  terras  of  the  inverted  functional  equation,  we  obtain 
F(x)-F(o)  =  J(VA1)fF[(x+l-a1)/a1]-F[(l-a1)/a1]) 

+  [ 1-^(A1+A2)  1  {F[(x-l+a2)/a2]-F[(a2-l)/a2]j 


[ft2-Al)/(2ai)  ]  -j^(x+l-a1)F[  (x+l-a1)/a1]-(l-a]L)F[  (1-^)/^] 

p  (x+1  -a,  )/a1 

F(u)  du  +K1+Kg 


-a. 


[(^2-^l)/(2a2)  ]^J-(x-1-hx2)F[  (x-l+a2  |/a2]-(l-a2)F[  (a2-l)/a2] 

(x-l4a2)/a2 


+a. 


’F(u)  du+K^+K^ 


F(x)-F(0)  =i[A1+^2+(A2-A1)(x+l-a1)/a1]F[ (x+l-a^/a^] 

+  [l-^(A1+A2)-i(^2-?v1)(x-l-Kt2)(x-lm2)/a2]F[(x-l+a2)/a2] 

r(x+l-a,  )/an  r  (x-l-HXp)/^ 

'  i(VV  /  F^du+^Ag-Aj.)  /  *  ^F(v)  dv 

J  _00  J  -00 

-  i[(A1+A2)-(A2-A1)(l-a1)/a1]F[(l-a1)/a1] 

-  [1-|(A1+A2)  +  i(A2-A1)(l-a2)/a2]F[(a2-l)/a2] 

+  |(A2-A1)(K1+K2)/a1  +  i(A2-A1)(K3+K4)/a2  . 


Imposing  the  condition  that  F(x)  =0,  x  <  -1  gives 

0  =  F(0)  -  i[A1+A2+(A2-A1)(l-a1)/a1]F[  (l-a^/c^] 
-  [1-|(A1+A2)  +  i(A2-A1)(l-a2)/a2]F[(a2-l)/a2] 
+  |(A2-A1)(K1+K2)/a1  +  i(A2-A1)(K5^-K4)/a2  . 


■  ;  "  ':o  -  \  :  .  -  O.  ■  h  \  /  ,  •  '  T\~'. 


: 


-  '  '  • 

et*.  • 
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Thus  the  limiting  distribution  functions  corresponding  to  the 
characteristic  functional  equation  (3.13)  satisfy  both  the  following  functional 
integral  equation 


(3.17)  F(x)  =  [^1+(^2"^l) (x+l)/(2a1) ]F[ (x+l-a1)/a1] 


+  [l-^2-(x-l)(A2-?v1)/(2a2)]F[(x-l-Ki2)/a2] 
p(x+l-a  )/a 

-i(W/,  F(u)  du  ,  0v,>\) 

J  (x-l+a „)/a 


and  the  auxiliary  conditions 


F(x)  =  0  ,  x  <  -1  ;  F(x)  =  1  ,  x  >  1  ; 

0  <  ,  a2  <  1  ;  0  <  <  A2  <  1 

Exactly  the  same  procedure  is  used  to  invert  the  characteristic 
functional  equation  (3.14)  to  determine  the  equation  satisfied  by  the 
corresponding  distribution  functions.  Under  the  inversion,  the  left 
member  becomes 


1_ 

2n 


(tp(t)[l 


e'1Xt]  /(it))  dt  =  F(x)  -  F(0)  . 


The  terms  of  the  right  side  of  (3.14)  become 

^  f°°  it(l-a-) 

i)  ^(A1+A2)  /  (Xc^t)  e  [1-e  ]/(it)}  dt 

J  —CO 


ii) 


=  J(A1+X2)[F[(x-l-Hx1)/a1]  -  F[(a1-l)/a1]]  . 

1  r00  -it(l-aQ) 

[l-i(A,+A  ))i-  /  (tpfagt)  e  -  (l-e‘1X  )/(it))  dt 

j  -00 


. 

\.( •  ! '  '  - 1  •  (  )•••■ 


x  <  0  -«  ;  >J 

■ 

•:C  :■ 


■ 
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=  [  [F[  (x+i-a2)/a2]  -  F[(l-a2)/a2J]  . 

-i(A  -A  )  r  °o  r  it  (l-a. )  .  h 

iii)  ~2^ -  2n  /_M  MalC>)e  [1-e'  X  ]/(it)|  dt 


Integrating  by  parts  gives 

-i(VV  r  i 


t)(« 


it(l-a1)  -it (x-l+a. ) 


-e 


1  /  c 

1  )/(it)] 


T=  -  oo 


it(l-a^)  -it (x-l+a^) 


-  ^  ^(ajt)  -^[(l-a^e  x  +(x-l+a1)e 


]/t 


[e 


it(l-a1)  -it (x-l+a^)  2 


A ,  -A  r  1-a, 


]/(it2)|  dt  |  i 

/'1-,'o  f  r1"  r  itCl-a.)  -it(x-l+a,)  "1 

--sqH--sr  LW)[e  -e  ]/(it)ldt 

.  roo  -it(x-l+a1)  -it(l-a1)  -it (x-l+a- ) 

+  ~  J  <P(alt)|(-xe  1  )/(it)  -  (e  1  -e  1  )/t 


where 


V*2  r  1 

— •J-(l-a1)[F[(x-l-Ki1)/a1]-F[(a1-l)/a1]]  +  h(x)|  « 


-it(x-l-Kx^) 


it(l-a1)  -it(x-l4a^)  • 


h(x)  =  ~  j'  cp(a;Lt)-j^-[xe  1  ]/(it)-[e  "-e  "  ]/t‘ 

-  f  ’it(X-1^1)  -itCx-l^)  .  •'U(X-1^1)  , 

h'(x)=!^  j  cp(<xxt )-J -  - - - -  +  xe  -  — - r -  rdt 


r00 

=  ~  /  9(axt)  e 

■J  _ on 


it 

-it(x-l+a^) 


dt 


(x/a^)  f  [  (x-l-KXj^  )/oc,1  ]  9 


provided  f  exists. 


■ 


* 


.  •  ■-  -  ;>-  !  , 


'  ■ 


....  •  • 


=  >  j 


I 
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Therefore 


1  fX 

h(x)  =  —  /  ufUu-l-w^/a  ]  du  +  ,  where  is  a  constant 


of  integration. 


h  (x)  =  (l-a^)F[  (x-l+a^)/a^  ]  + 


(x-l+a..  )/a. 

v  f(v)  dv  +  f 


=  (1~cc^)F[  (x-l+a^)/a^]  +  (x-l+a^)F[  (x-l+a^ )/a^ ] 


-a. 


(x-l+a^ )/a^ 


F(v)  dv]  +  . 


Thus  the  inverted  form  of  iii)  is 


A  -A 


-  X(x-l+a1)F[  (x-l+a1)/a][]  +  (l-^ )F[  (a^l  )/a1  ]  -  f  ( 
]_  L  J  .00 

F(u))du  +  Kx  +  Kgj-  . 


(x-l+a^) 


i(A--Ap)  .  o~  r j 


[  1-e”  ]/(it)j-  dt 


Again  integrating  by  parts  we  obtain 


i(Ai-A2)  r  1  ,  ^  -it(l-a2)  -it(x+l-a2) 


2a, 


{  51  [<P(V)(' 


-e 


)/(it)] 


t=  -00 


,  O00  r  -it(l-a  ) 

+  —  /  cp(a0t)4  [  (l-a2)e  -  (x+l-a2)e 


-it(x+l-a2) 


2jt  J  ^  T'  2 
y  -00 


]/t 


-it  (l  -aQ)  -it(x+l-a0) 

+[e  2  -e  2  ]/(it£) 


■ 


•  ■ 


o 
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ArA?  r  1_ap  r00  r  -it(i-a  )  -it(x+i-a  ) 

|cp(apt)[e  ^  -e  ]/(it)f"dt 


2a2  l  2*  J-oo  1T'~2' 


+ 


oo  r 


-it (x+l-a  )  -it(l-ap) 

xe  2  )/(it)+(e 


-it (x+l-a  )  0. 

■  -  2  )A2 


^|-(l-ag)  jj[  (x+l-ag)/a2]  -  F  [  ( 1  -ag  )/ag  ]  j-  +  g(x)  , 


where  g(x)  =  g-  /  q>(a,2t)[(xe 

w/  -00 


-it(x+l-a2)  -it(l-a2)  -it(x+l-a0) 


)/(it)  +(e  -e 


2  )/t2]  dt 


s’W  =  a?/ 

J  -00 


e-it(x+l-“2)  -it (x+l-a  )  e-iAK+l-a2) 
— -  -  xe  -  - 


it 


it 


dt 


2*  /  (Ka2t)e 

-CO 


■it (x+l-a  ) 


dt 


=  -(x/a?)  f  [  (x+l-ag)^^]  ,  (if  the  derivative  f(x) 


of  F(x)  exists ) . 


Thus 


i 

g(x)  =  -  —  /  uf  [  (u+l-a2)/a-2]  du  +  K  ,(K  is  a  constant  of 

a2  J  -00  ^  ^ 

integration) , 


( g (x ) )  =  (l-a2)F[(x+l-a2)/a2]  -  ag  J' 


(x+l-a2)/a2 


uF(u)  du  +  K  , 


=  (l-a2)F[ (x+l-a2)/a2]  -  (x+l-a2)F[ (x+l-a2)/a2] 

^(x+l-a2)/ap 


+  a. 


F(u)  du  +  . 


-00 


'  (  ):■') 


Thus  the  inverted  term  is 


^1 f 

2^ — -j  -  (l-a2)  [F[  (x+l-a2)/a2]-F[  (l-a2)/a2]  ]+(l-a2)F[  (x+l-a2)/a2J 

r  (x+l-a  )/a? 

-  (x+l-a2)F[  (x+l-a2)/a2]-Kx2  /  F(u)  du  +  K  +K^ 


2a^ 


■|(l-a-2)F[(l-a2)/a2]  -  (x+l-a2)F[  (x+l-a2)/a2] 

r  (x+1-a2)/a2 


+  a. 


-00 


F(u )  du  +  K  +  K1( 


Collecting  the  terms  of  the  inverted  equation  we  obtain 


F(x)  -  F(0)  =  {F[(x-l4a1)/a1]  -  F[  (c^-l)/^ ] ) 

+  [l-i(A1+?v2)  ]  {F[(x+l-a2)/a2]  -  F[  (l-a2)/a2 ] } 

+  — -|^(x-1-kx1  )F[  (x-l+a1  )/oc-1  ]  +  (l-a1)F[  (a1-l)/a1l 


-  a. 


(x-l+a^ )/a^ 


F(u)  du  +  +  K, 


2ag  I  '  2' 


(l-aQ)F[ (l-a2)/a2]  -  (x+1 -ap)F[ (x+l-ap)/ap ] 


2"  2 


p (x-l+a )/a 

+  a2  /  ^  ^  F(u)  du  +  K  + 

J  -00  ^ 


Imposing  the  condition  F(x)  =0  ,  x  <  -1  gives 


F(0)  =  i[A1+A2  -  (A1->v2)(l-a1)/a1]F[  (a1-l)/a1] 

-  U-i(\+A2)  "  i(^1-\)(1“a2)A2]F[(1-a2Va'2^ 

+  i(A1-A2)[(K1+K2)/a1  +  (K5+K4)/a2]  . 


Thus  for  >  ^2  the  limiting  distributions  satisfy  the  functional  equation 


.. 
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(3.18) 


F(x)  =  [A1+|(\-A2)  (x-1  )/ax  ]F[  (x-1-kx^  )/^  ] 


+  [l-^2"4(^i“^2^X+1^a2^F^  (x+1_a2Va2] 

r(x+l-a  )/aP 

+  J  F ( u ) du ,  (\  >  >2,  \  /  1,  /  0)  , 

(x-l+a^  )/a. 


and  the  auxiliary  conditions 


F(x)  =0,  x  <  -1  ;  F(x)  =  1  ,  x  >  1  ; 


0  <  al>  a2  <  1 


;  o  <  ~h2  <  ~h1  <  l  . 


§  3.1+  Additional  Properties  of  the  Limiting  Distributions 

The  characteristic  functions  of  the  untransformed  limiting 

distribution  functions  satisfy  the  equation 

iA  (i-a  )0  iA  (1-a  )0 

cp(e)  =  e  cp(a20)-(i/a1)e  ^{9(0^0)] 

i^  (1-a  )0 

+  (i/a2)e  2q  [9(^2®) 1 • 

This  equation  may  be  re-written  in  the  form 

iA  (l-co  )0  iA  (1-a  )0  iA  (l-a  )0 

(3.19)  9(0)  =  e  ^  ^  cp(a20)-ie  cp'fa^J+ie  ^(a^), 


wh 


ere  t>”K9>  - 


r^2e)  =  WV” 


and 


Since  the  characterstic  functions  have  the  property  that  cp(0)  =  1 


V„  -  d/in)  =0 

dd 


tin 

where  V  is  the  n  moment  of  the  desired  distribution  function,  we 
n 
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seek  to  derive  the  moments  of  the  limiting  distribution  from  the  characteristic 
equation  (3.19)»  For  (9=0,  equation  (3. 19)  reduces  to 

0  =  (-i  +  i)  cp'  (0)  . 

Thus  <p'  (0)  is  indeterminate. 

Differentiating  (3. 19)  with  respect  to  6  gives 

(1-a  )(9  i\  (1-a  )(9 

q)'(@)  =  iA2(l-a2)e  9(cc20)  +  [a2»^2(l-a2)  ]e  cp'^e) 

iA  (l-a..)0  iA  (l-a  )(9 

+  A^(l-a^)e  cp' (a^0)-ia^e  ~  cp"  (a^S  ) 

iA  (1-a  )(9 

+  ia^e  9"  (a20)  . 

For  £9  =  0  the  above  equation  becomes 

-iA2(l-a2)  =  [a2  -  +  ^1(l“a1)  "  1 (O)  +  i(a2-a1)cP"  (0)  . 

This  equation  cannot  be  solved  in  general  since  there  are  two  unknowns  cp 8  (0 ) 
and  cp"(o)  .  However,  if  a^  =  a2  >  t^ie  coefficient  of  cp"  (0)  vanishes. 

Thus  replacing  and  a2  by  a  ,  we  have 

CP '  (0)  =  iAg/(l  +  -  \)  ,  (a1  =  a2  =  a  )  . 

Lemma  3 . 1  The  functional  differential  equation  satisfied  by  the 

characteristic  functions  of  the  limiting  distribution  of  response  probabilities 
specifies  unique  moments  for  =  a2  =  a. 

Proof  Imposing  the  equal  alpha  condition  on  (3-19)  gives 

iA  (l-a)(9  iA  (l-a)£9  iA  (l-a)£9 

(3.20)  cp(@)  =  e  ^  cp(a0)  +  i[e  ^  -  e  ]cp»(a<9)  . 


i-  '  0 


v  '•  -  ■  p 
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Differentiating  (3.20)  n  times  with  respect  to  Q  then  imposing  6=0 

/n . 1 \  (n) 

causes  the  coefficient  of  cp'  '(O)  to  vanish.  This  leaves  cp'  '(0)  as 

fk) 

a  function  of  cp'  '  (0)  (k=0, 1,  .  . .  ,n-l) .  Thus  the  moments  {V^}  maY 

calculated  in  order  of  increasing  n. 


Direct  verification  shows  that  the  result  of  lemma  3 »1  also  applies 
to  the  transformed  characteristic  equations  (3.13)  and  (3*14). 

It  has  been  shown  above  that  the  limiting  distributions  of  the 
transformed  response  probabilities  of  the  Subject-Controlled  Model  satisfy 
either 


(3.21) 


F(x)  =  [A1  +  (x+l)(^2-7v1)/(2a1)]F[(x+l-a1)/a1] 


+  [1  -  +  (l-x)(A2-7\1)/(2a2)]F[(x-l4a2)/a2] 

r(x+l-a  )/a 

-  |(*2“\)  J  F(u)  du  ,  (A2  >  Ax)  ; 

(x-l+a2)/a2 


or 


(3.22) 


F(x)  =  [A1  +  (x-l)(^1-7v2)/(2a1)]F[(x-l4a1)/a1] 


+  [l-^2  -  (x+l)(^1-^2)/(2a2)]F[(x+l-a2)/a2l 
n (x+l-a  )/a 

+  )  J  F(u)  du  ,  (^  >  ^  l> 

(x-l+x1)/a1  ^  Q) 


with  the  auxiliary  conditions 


(3-24) 


F(x)  =  0  ,  x  <  -1  ;  F(x)  =  1  ,  x  >  1  ; 
0  <  al  ,  a2  <  1  ;  0  <  \  -  1 


. 
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Theorem  5.1  The  limiting  distribution  functions  satisfying  the  functional 
systems  (3 .21 ) , ( 3.23) , (3 -24)  and  (3.22) , (3.23) , (3*24)  are  unique  if  a^=a^=a0 

Proof .  Lemma  3*1  proves  that  the  characteristic  functions  are  unique 

for  ai=a2=:a"  Thus  by  the  uniqueness  theorem  of  characteristic  functions, 
there  can  be  only  one  distribution  function  that  satisfies  the  above  functional 
systems  (c.f.  Theorem  2.3). 


Theorem  3.2  For  a^+a^Cl  ,  the  distribution  functions  satisfying  the  above 
functional  systems  are  purely  singular  functions  which  are  constant  everywhere 
except  perhaps  on  a  non-denumerable  set  of  points  of  Lebesgue  measure  zero. 


Proof .  The  arguments  of  (3.21)  and  (3.22)  are  of  the  form 

(x+l-a')/a'  ,  x,  (x-l+a*)/a* 

(For  (3.21)  read  for  a'  and  a 2  for  a*,  for  (3.22)  read  co2  for  a? 

and  for  a*).  In  the  domain  (-1,1)  the  arguments  are  ordered  from  largest 

to  smallest  as  written  above.  For  a^+a^l  >  we  know  from  theorem  2.1  that  there 

exists  an  interval  I1  1=[2a' -1, l-2a*]  such  that  (x+1 -a" )/cl1  >  1  and 

i ,  i 

(x-l+a*)/a*  <  -1  simultaneously. 


Consider  (3.21)  first.  Letting  xel^  ^=[2a^-l, l-2a2]  and  applying 
the  boundary  conditions  (3.23)  gives 

(x+1 

”al  )/al 


r 

F(x)=>v1+(x+l)(A2-A1)/(2a1)-^(A2-A1)  / 

-1 


Therefore 


=A1+(x+l)(A2-A1)/(2a1)-i(A2-A1)|~  J F(u)du+  ^ 

F(x)=A2-J(A2-A1)  J  F(u)  du  ,  (x  e  I1^1) 


F(u)  du  , 
(x+1  )/a. 


du 


The 


above  equation  shows  that  F(x)  is  constant  for  all  x  el^  Let  this  constant 


be  F^  ^  so  that 


(3.25)  Fx  /  F(u)  du  • 

i  J  -1 

Since  the  value  of  J'  F(u)  du  is  unknown,  (3.25)  does  not  yield 


. 

. 

«  '  "  -  -  ‘^e"  .2.3)  smg;)  v; 


«. .)  &n.  •  ~.2.£)  So  aitnami^oo  s>1T 
■\(*£  ■  ,  -\{  !v>Xj  -<} 
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the  explicit  value  of  F.  in  general. 

1,1 

Using  the  notation  introduced  in  Chapter  II  (p.  2k  f.f),let  G1 

denote  the  gap  (-1,  2tt^-l)  ,  and  ^  denote  the  gap  (l  -  2^,1). 

From  the  interval  I,  ,  ,  two  new  disloint  intervals  where  the  distribution 

1,1* 

function  takes  a  constant  value  may  be  found  by  the  operators  (2. 37) 

and  The  value  of  the  distribution  function  on  the  interval 

I2  1  =  [a1(2a1-l)-(l-a1) ,  (l-2a2) - ( 1-a^ ) ]  is  given  by 

-  p  2a1  - 1 

F(x)  =  [A1+(x+l)(>2-A1)/(2a1)]Fl^1  -  J  F(u)  du 

p  ( x+1  -a^ ) /a^ 


2a^-l 


(  X  e  I2A)  . 


Fl,l  du 


Simplifying,  we  have 

F(x)  =  [A1  +  a1(^2“Al)]Fljl  "  J  L  F(u)  du  * 


, 2a^-l 


Since  the  solution  on  ^  is  also  a  constant,  let  ^2  1  =  €  F2  1^ 

so  that 

p2a.  -1 

0.26)  f2>1  =  [Xj  +  ■  ?(w  f(u)  du  • 


Similarly,  the  distribution  function  on  the  interval 


I2  2  =  [a2^ai”1 )  +  a2^1"2a2^  +  ^1_a2^ 

is  given  by  F(x)  =  f2  2  ,  (x  eI2,2^  where 

(3.27)  F2)2  =  +  [l-\,  +  a2(V\)]F1(1  -  *(VV  J  F(“)  du 


.  ,  '  :  '  ■;  "■  ;  i  j  a?  •  r 


ft  J  1v 
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In  this  problem  a  set  of  disjoint  intervals  over  which  the 
distribution  function  is  constant  is  generated  in  exactly  the  same  manner 
as  the  intervals  of  the  Experimenter  Model.  Thus  the  intervals  are  disjoint 

and  may  be  generated  in  a  series  of  stages  such  that  at  the  ntk  stage  2°  * 

n_2 

new  intervals  are  generated  from  the  2  intervals  of  the  n-1  stage  by 

the  transformations  and  B^.  We  now  proceed  by  induction  to  show  that 

the  value  of  the  distribution  function  in  the  interval  1^  ^  is  given  by 

L  .  . 
r  n-l,i 


(5.28)  Fn>.  =  [\4(VV(1+Ln-1(i>lFn-l,i-i{VV  / 


F(u)  du 


n-2> 


(n=2,3,4, ...,  i  =  1,2,. ..,2  "  ) 


or  by 

(3.29) 


Fn,i  ■  *2  +  [!-Vi(A2-  V(U  .  .  9n-2-DF  ,  .  --2 

n-1, i-2  n-1, i-2 


,  •  on-2 

n-1, i-2 


n— 2 

F(u)  du  ,  (n=2, 3, . . . > i=2  +1, 

2n-S+2,...,2n-1)  , 


where  [L  ,  U  ]  is  the  interval  I  and  F  is  the  value 

n-l,i  n-l,i  11  •L*-L  ’ 


>f  the  distribution  function  over  1  .  * 


n,  l 


The  induction  hypothesis  has  been  proved  for  n=2  .  Suppose 
that  it  holds  for  n  =  3,4,...,  k-1  ,  For  n=k,  the  intervals 
I  . (i=l,2, . . .,2k'2)  are  generated  from  I  .  by  the  transformation 

k,  i  ’ 

Bl.  Thus  x  e  I.  .  for  all  (x+l-a,)/^  e  Ik-1  r  Now  let  the  value 
of  the  distribution  function  over  Ik_1^i  be  Fk_^i‘  From  equation 
(3.21)  and  the  lower  boundary  condition  (3*23)  we  §et 

r (x+l-a] )/aj 

F(x)  =  [A1  +  (x+l)(V\)/(2a1)]Fk.i,i-i^2-\)  J  F(u)  du, 

-1 


(X  £  Ik,i)  ’ 
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Lk-l,i  /x+1-ai)/“i 


-1 

L 


F,  du 

k-1 ,  i 


k-1 ,  i 


-  p  k-1 , i 

=  [A1+(x+l)(>v2^1)/(2a1)]Fk_1  .^(^g-^)  /  F(u)  du 

’  u  -1 


Jk-l,i 


=  [\-<i(VAi)(l4lfc.i,ii  Fk-i,ri(VV/_1  F<u)  du  >  (xeIk;1)- 


Since  F(x)  is  a  constant  over  I.  ,  it  may  be  denoted  F.  .  Thus 

K  y  1  K  y  J_ 


Fk,i  =  [\4(V\}(1+Lk-l,i)lFk-i,i ’  F<u)  du  > 


(x  e  Vi5 


This  completes  the  proof  of  (3.28). 


The  same  type  of  argument  may  be  used  to  prove  equation  (3.29). 

In  this  case,  the  intervals  1^  ^  (i=2k  2+l,  2k  2+2,  ^)  are  generated 

from  I  .by  the  transformation  B  .  Thus  for  x  e  I,  .  , 

k-l,i-2k-2  k_2  2 

(x-l-toO/a^  el  (i>2  ).  For  simplicity  of  notation,  let 

d  k-l,i-2 

k-2 

i'  =  i-2  ).  From  equation  (3.21)  ,  the  upper  boundary  condition  (3.23) 

and  the  induction  hypothesis  we  have 


F(x)  =  A1+(x+l)(A2-A1)/(2a1)+[l-A2+(l-x)(A2-A1)/(2a2)]Fk_1^i, 

-  rUk-l,i'  pi  p  ( x+1  -a.  )/a, 

-  i(VV  /  Fk-i,i'du  +  /  F<u)du  +  / 

(x-Ua  )/a  Uk-1, i'  1 


du 


■ 
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=  ^1+(x+l)(A2-X1)/(2a1)+[l-A2+(l-x)(A2-A1)/(2a2) ]Fk_1  £| 

‘  ^2'\^  (Uk»i^  i » "1“(x“1)/Q'g)Fk_i^  i  t  +  ^  F(u )du  +  (■x.+l)/a,l-2. 


U 


k-1,  i 


Therefore 


F(x)=>^+[1-A2-J(A2-^1 )(Uk_1>i,-l)]Fk_1^i,-|(A2-A1)  J  F(u)  du, 


U 


k-1, i’ 

/  T  •  ok -2  ,  ok-l\ 

(x  e  Ik  i,  1=2  +1, ... ,2  ) . 


Jc-2, 


Since  F(x)  is  constant  over  the  intervals  I,  .  (i  >  2  ) 

k,  i 

write  F(x)  as  F,  .  so  that 

k,  i 

Fk,i  ■  %  F(u)du 


k-1, i ' 


Thus  we  have  shown  that  the  limiting  distributions  which  satisfy 
(3.21),  (3.23),  (3.21+)  for  <Xl+a2  <  ^  3re  constant  otl  a  set  °f  intervals 


J-  •  • 

n,  i 

(3.22) 

(3.30) 


The  same  method  of  proof  may  be  applied  to  the  functional  system 
(3.24)  for  ai+a2  <  1  to  skow  that 

Fk,i -  /_x  ’  F(u)  du  ’ 


and 


k-2. 


(x  e  Ife  .,  i=l,2, . . .,2  ) 


(3.31)  Fk^  =  l-A1+[A14(A1-A2)(Uk.ljil-l)]Fk_lji,+l(A1-A2)  J  P(u)  du 

k-l,i’ 


(x  e  Ik  t,  i=2k“2+l,  2k'2+2, . . .,2k‘1) . 


**1 


1 
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In  theorem  2.1  it  was  shown  that  the  set  of  intervals  that  are  generated 
by  the  transformations  B^  and  B0  are  disjoint  and  cover  the  domain 
(-1,1)  except  for  a  non-denumerable  set  of  points  of  Lebesgue  measure 
zero . 


In  the  case  of  a^+a^  >  1  ,  almost  nothing  is  known  about  the 
limiting  distribution  functions.  The  trial  functions  which  satisfied  the 
functional  equation  and  boundary  conditions  in  the  Experimenter-Model  do 
not  satisfy  the  functional  integral  equations  of  this  model. 


f  :  --  .  •  •  ’  ' 


■ ;  - 

.  ;  ’• 
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CHAPTER  IV 

NUMERICAL  SOLUTIONS  FOR  THE  TRANSFORMED  LIMITING 
DISTRIBUTIONS  OF  THE  EXPERIMENTER-CONTROLLED  MODEL. 

§4.1  Introduction 

Numerical  approximations  of  the  limiting  distributions  of  response 
probabilities  are  of  practical  value  for  both  applied  statistics  and  pyschology. 
This  may  be  especially  true  in  the  purely  singular  cases  of  >  1  where 

numerical  tables  may  be  the  only  method  of  characterizing  the  distribution. 

In  this  chapter  a  numerical  method  is  described  for  tabulating  the  limiting 
distributions  by  a  numerical  solution  of  the  functional  system  (2.26)  -  (2.28). 
This  procedure  is  primarily  intended  for  the  cases  a^+a^  >  1  since  the 
distributions  for  a^-Kt ^  <  1  can  be  uniquely  determined  by  the  algorithm 
of  theorem  2.1.  Unfortunately  a  general  error  analysis  of  the  numerical 
procedure  is  impossible  at  this  stage  due  to  the  lack  of  detailed  information 
about  the  nature  of  the  limiting  distributions  for  o^-fct ^  >  1.  Thus  we  have 
no  general  proof  that  the  procedure  will  converge  to  the  true  distribution. 

The  best  that  we  can  do  is  to  compare  the  numerical  results  with  the  known 
properties  of  the  distribution  in  the  special  cases  where  this  can  be  done. 

An  extremely  favourable  comparison  is  obtained  between  the  numerical  results 
and  known  absolutely  continuous  solutions,  centro-symmetric  distributions 
and  a  previously  published  Monte-Carlo  approximation  given  in  [1].  Thus 
we  have  a  strong  indication  that  the  numerical  procedure  yields  accurate 


solutions . 
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§  4,2  Numerical  Procedure 

The  major  problem  is  to  devise  a  fast  and  accurate  method  for 
tabulating  the  numerical  values  of  the  limiting  distribution  at  a  set  of 
points  in  the  domain  (-1,  l).  Since  the  functional  system 

(k*l)  F(x)  =  J^Ft  (x+l-a^/c^)  +  rt^[(x  -l-KX2)/a2)  , 

(4.2)  F(x)  =  0  ,  x  <  1  ;  F(x)  =1,  x  >  1  * 

(4.3)  0  <  ctp  a2  <  1  ;  0  <  )  *2  -  1  >  «i  +  itg  =  1  , 

defines  a  unique  distribution  function  (Theorem  2.3),  we  may  reasonably 
expect  that  a  numerical  solution  is  possible. 

To  begin  the  numerical  procedure  we  first  define  a  set  of  tabular 
points  {x^}  (k=0,l,2, . . ,,2N)  lying  in  the  closed  region  [-1,  1],  For 

simplicity,  we  will  initially  choose  the  tabular  points  to  be  equi-spaced 
such  that  the  tabular  interval  h  =  ^-x^I/N.  On  this  tabular  set  {x^} 

the  functional  equation  becomes 

(4.4)  F(kh-l)  =  [Ich-^J/a^^  ]  +  jtgF  :  (kh-2-«x2)/a2  1  •  (k=0,l,2, . .  .,2N)  . 

In  general,  the  arguments  (kh-a^)/a^  and  (kh-2+a2)//a2  will  not 
lie  on  any  member  of  the  chosen  tabular  set.  Thus  (4.  4)  defines  a  system  of 
2N-1  equations  in  6N-3  unknown  values  of  F,  Imposing  the  boundary 
conditions  (4.2)  reduces  the  number  of  unknowns.  The  lower  condition  gives 

(4.5)  F(kh-l)  =  rt^F  j(kh-a1)/<x1  1  ,  (k*l,2,3, . . . ,m)  ; 


... 


i 

. is.  .  f  • ;  .cl  rr,r 

■ 

' 

■ 

■ 
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where  mh-1  <  l~2a  and  (m+l)h-l  >  X-2a,, 

d  d 

The  upper  boundary  condition  implies 

(4.6)  F(kh-l)  =  3^  +  itgF  ;(kh-2-Kt2)/a2  )  ,  (k=n,n+l, . . .  ,2N-1 )  ; 

where  nh-1  >  2a ^  - 1  and  (n-l)h-l  <  2a.. -I  . 

Because  there  are  no  more  linearly  independent  conditions  involving 
only  the  6N-3-m-n  values  of  F  already  given  in  (4.4)  -  (4.6)  ,  the  only 
feasible  approach  is  to  interpolate  between  the  tabular  values  to  approximate 
the  off-tabular  points.  In  this  problem  only  the  function  values  at  the 
tabular  points  are  at  our  disposal.  This  makes  it  rather  inconvenient  to 
apply  the  most  accurate  interpolation  methods  like  the  Bessel  and  Everett 
Interpolation  formulas,  since  these  procedures  are  expressed  in  terms  of 
the  differences  of  the  function  values.  Furthermore  the  interpolating 
factor  p  (x-ph  =  x^)  will  probably  be  different  for  each  point  to  be 
interpolated.  This  would  necessitate  calculating  a  whole  set  of  Bessel  or 
Everett  coefficients  for  every  off-tabular  point. 

Thus  we  try  the  simpler  but  less  accurate  Lagrange  Interpolation 

formula. 

Having  chosen  the  interpolation  formula,  the  next  problem  is  to 
determine  the  order  of  interpolation  that  will  give  the  best  results.  Be¬ 
cause  many  of  the  solutions  are  known  to  be  purely  singular,  a  high  order 
formula  may  introduce  instability.  Even  the  known  absolutely  continuous 
solutions  have  discontinuities  in  their  higher  derivatives.  In  the  latter 
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case,  a  high  order  interpolation  formula  would  either  accentuate  the  effect 
of  the  discontinuities  or  smooth  them  over  and  hide  them.  For  these  reasons 
a  four  point  Lagrange  interpolation  formula  was  tried,  with  the  initial 
scheme  to  choose  two  tabular  points  on  each  side  of  the  off-tabular  argument. 


After  a  number  of  numerical  experiments,  the  following  scheme 
proved  to  be  far  superior  to  any  of  the  other  variants  tried.  Given  an  off- 
tabular  point  x'  lying  between  x^  and  x^+^  (k  ^  0,1,  2N-2,  2N-l)  F(x') 
was  approximated  by  the  formula 


(4.7) 


F(x’)  = 


■x  )(x'-x  ) 

k+1  k-t  2 


(x’-XfcH: 

(\-r\)(\-r\+i)(xk-rxk+2) 


F(xk-i) 


(x '-x,  )(x'-x  )(x'-x  ) 

.  k-1  k+1  k±2  „/  \ 

(Vxk-i)(V\+i)(V\+2) 

^Xk+l'Xk-l^\+l"Xk^Xk+l"Xk+2^  k+2 


For 


-1  <  x'  <  x  the  interpolation  formula  was  changed  to: 

d. 


(x’-x2)(x'-x^)(x'-x^)  ^  ^  (x,-x1)(x'-x5)(x’-xl|_) 

(x1-x2)(x1-x5)(x1-x^y  F^X1^  +  (x2“xl )  (x2-x3)Tx^)  F^X2') 


(x'-x1)(x,-x2)(x'-x1^)  (x'-x1)(x'-x2)(x,-x  ) 

(x3-x1)(x3-x2)(x5-xu)  F(x3)  +  TvxiKx4"X2^Xl+“Xy  F^Xl|) 


(4.8)  F(x') 


' 


' 


T*  ”  •  •  ..  -4 
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And  for 

(4.9) 


the  form 


(4.10) 


x__T  ^  <  x'  <  1  ,  the  formula  used  was: 

2N-2  * 

w  ,x  _  (X,-X2N-2)(X'-X2N-1)(X,-X2N)  w  * 

^  X2N - 5 " X2N -2  ^  X2N - 3 " X2N - 1 ^  X2N -  3  ~ X2N  ^  2N"5 


(X,-X2N;3)(X,~X2N-1KX,~X2N)  p(x  } 

^X2N-2'X2N-3^X2N-2"  2N-1^X2N-2'X2N^  2N~2 


^X,-X2N-3)(X,~X2N-2KX,-X2N}  r(::  } 

^X2N-l"X2N  -  3  ^  X2N  - 1 "  X2N  -  2  ^  X2N  - 1  “  X2N^  2N~1 


^X,~X2N-3|l/x'“X2N-2^X,~X2N-l^  p,  x 
^  X2N  '  X2N  -  3  ^  X2N  "  X2N  -  2  ^  X2N  “  X2N -1^  2N 


To  set  up  the  system  of  linear  equations  in  F(x^) (k=l,2, . . .  2N-1 ) 
of  (4.4)-(4.6)  was  changed  to 


■F(kh-l)  +  F  |[(kh-a^)/a^  =  0 


(k— 1  , 2? , . , ?m)  } 


-F(kh-l)+jt^F  [(kh-a^/a^  j+^rt^F  [ ( kb  - 2 ) /cc^  ]  =  0 

(k=m-hl,m+2,  .  . .  ,n-l ) 

-F(kh-l)  +  jtgF  [  (tch-2-xx2)/a2  j  =  ;  (k=n, n+1, . . . ,2N-l) 


The  final  system  of  equations  has  the  following  form. 
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where 


Cij 


are  the  Lagrangian  Interpolation  coefficients  for  the  argument 


(lch  -  a1)/a1  , 


and  c .  . 
ij 


are  the  Lagrangian  Interpolation  coefficient  for  the  argument 


(kh-2-KX2)/a2  . 


A  Fortran  program  was  written  to  set  up  and  solve  the  above 
system  of  equations  given  any  set  of  a^,  a2,  and  Ttg,  along  with  a 

tabular  interval  h.  The  final  matrix  equation  was  solved  using  a  Gaussian- 
elimination  procedure  which  chooses  the  best  pivots  of  each  column. 


The  exact  procedure  used  to  set  up  the  system  of  equations  for 
tabulating  F(x^)  (k=l,2, . . . ,2N-l)  is  best  explained  by  the  following 
ALGOL  program. 

procedure  Experimenter  Limiting  Distributions; 

real  array  A[l:  M,  1:  M+l],  x[0:  M+l],  F[0:  M],  c[ljl+]  ; 
real  h,  alfa  1,  alfa  2,  pi  1,  arg  1,  arg  2,  swtch,  pi  2,  pi  ; 
Boolean  leap,  jump,-  Integer  i,  j,  k,  max,  lead; 

comment  This  procedure  computes  the  values  of  the  transformed 
limiting  distribution  function  of  the  Experimenter-Controlled  Model 
at  a  set  of  points  x[k]  in  the  domain  (-1,1).  These  values  are 
determined  by  a  system  of  linear  equations  in  the  unknowns  F(x[k]) 

procedure  Lagrange  (lead,  max,  arg,  x,  c)  ; 

real  d  ;  integer  i,  j,  k,  £  ; 


value  lead,  max,  arg  ; 
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comment  Lagrange  is  a  sub -procedure  which  calculates  the  Lagrange 
Interpolation  coefficients  c[i]  given  the  value  of  the  abscissa 
' arg',  the  set  of  tabular  points  x[i]  ,  and  the  index  'lead' 
which  defines  the  set  of  points  to  be  used  in  the  interpolation. 
The  value  'max'  determines  if  the  upper  boundary  is  used  in  the 
interpolation. 

begin 

for  i:  =  1  step  1  until  4  do 
begin 

d:  =  1.0  ;  k:  =  lead  +  i  ; 

for  j:  =  1  step  1  until  4  do 

begin 

if  i  =  j  then  go  to  cycle  else 
begin 

i:  =  lead  +  j  ; 

d:  =  d  x  (erg  -  x[ &] )/(x[k] -x[ &]  ] 

end 

cycle:  end  loop  with  index  j  ; 
c[i]  ,  =  d  ; 

end  main  loop  for  calculating  c [ i ]  \ 
if  lead  +  3  =  max  then  c[4]  :  =  -  c|.4] 
end  Lagrange  This  completes  the  sub-procedure  Lagrange} 
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begin  Comment  The  following  program  segment  is  the  main  body  of  the  procedure 
which  sets  up  the  matrix  equation; 
start i  read  (h,  alfa  1,  alfa  2,  pi  1  )  ; 

max:  =2.0/h  ;  pi  2:  =  1,0  -  pi  1  ; 

for  i:  =0  step  1  until  max  +  1  do  x[i]:=hxi— 1.0; 

arg  1:  =  1.0  -  alfa  1;  arg  2:  =  alfa  2  -  1.0;  swtch:  =  1.0  - 

alfa  1  -  alfa  2  ;  leap:  =  true  ; 

for  i:  =1  step  1  until  max  -1  do 

begin  comment  This  loop  determines  the  i*^  row  of  the  matrix; 
for  j:  =  1  step  1  until  max  do  A[i,  j]:=0.0; 

A[i,i]:  =  -1.0  ;  pi:  =  pi  1  ;  jump:  =  true  ; 
if  leap  them  go  to  execute  else  go  to  skip  ; 
execute:  arg*  =  (x[i]  +  arg  l)/alfa  1  ; 

if  arg  >1.0  then  go  to  set  else 
begin 

search  :  for  j:  =  1,  j+1  while  x[j]  <  arg  do  \ 

if  x[ j ]  =  arg  then  A[i,j]:  =  pi  +  A[i,j] 


else 

if 

j 

< 

2  then  lead :  -  0 

else 

if 

j 

> 

max  -  2  then  lead :  =  max-4 

else 

lead : 

= 

J  -  3  ; 

end 

Lagrange  (lead,  max,  arg,  x,  c)  ; 

for  k:  =  1  +  lead  step  1  until  4  +  lead  do 


A[ i,k] :  =  A[i,k]  +  pi  x  c[k-lead]  ; 


. 


-  .. 

•  1  ! : 


' 

■ 
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go  to  reset; 
leap:  =  false ; 

A[i,  max]:  =  -pi  > 

if  jump  then  begin  if  x[i]  <  swtch  then  go  to  restart  end 
else  begin  comment  this  part  re -initializes  to  interpolate 
for  the  smallest  argument; 
pi:  =  pi  2  ;  arg:  =  (x[i]  +  arg  2)/alfa  2  ; 

jump:  =  false ; 

go  to  search; 

end 

restart:  end  This  ends  the  main  loop  for  setting  up  the  matrix  equation; 

GAUSEL  (max  -1,  A,  F)  ; 

Comment  GAUSEL  is  a  Gaussian  elimination  procedure  for  solving  sets 

of  linear  equations.  Given  matrix  A  of  order  (max  -l)X(max), 
this  procedure  solves  the  set  of  equations  by  Gaussian  Elimination 
with  best  pivots  and  leaves  the  solution  in  the  column  vector  F; 
F[0] :  =  0.0  ;  F(max):  =  1.00; 

Punch  (alfa  1,  pi  1,  alfa  2); 

Punch  (x[i],  F[i],  i:  =  0,  max); 

Comment  The  above  segment  merely  punches  out  the  values  of  the  parameters 
then  a  table  of  x[i]  and  F[i]; 
go  to  start; 


set : 
skip : 
reset : 


end  Experimenter  Limiting  Distributions; 


’ 
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§  4.3  Accuracy  of  the  Numerical  Results. 

The  Fortran  program  for  executing  the  algorithm  of  §4.2  was  run 
on  the  I.B.M.  1620  system  of  the  University  of  Alberta.  The  calculations 
were  done  using  hardware  floating  point  carrying  an  eight  digit  mantissa 
and  two  digit  exponent.  A  number  of  test  cases  were  calculated  to  compare 
the  numerical  answers  with  values  which  could  be  determined  from  theoretical 
considerations.  All  the  test  cases  indicated  a  very  satisfactory  accuracy- 
more  than  adequate  for  any  practical  problem.  However,  because  of  the  nature 
of  this  problem,  no  error  bounds  on  the  numerical  process  itself  have  been 
determined . 


The  simplest  of  the  test  cases  depended  upon  the  property  of  centro- 
symmetry  of  the  transformed  limiting  distributions  for  a]_  =  a2  =  a  > 


rt,  =  = 


_  l 


.  This  property  may  be  proved  from  a  result  of  McGregor  and 


Hui  ([2]p.l282).  These  authors  showed  that  in  the  equal  alpha  case,  the 
characteristic  function  of  the  transformed  limiting  distribution  reduced  to 

oo  -i@  ( l-a.  )a, n  .  \  n 

(v.ii)  <p(e)  =  n  tv  «2eie(1-a)a  ]  . 

n=o 


Choosing  jt^  =  n  =  T-  ,  the  above  characteristic  function  may  be  written  as 
an  infinite  product  of  cosines.  Thus  the  frequency  function  must  be  symmetric 
about  zero  and  the  distribution  function  centro-symmetric  about  (0,i) .  The 
first  test  consisted  of  tabulating  the  distribution  functions  corresponding 
to  =  rt2  =  a  taking  the  values  .55(05)  ’95*  The  numerical 

solutions  for  the  tabular  interval  h=0.1  are  displayed  in  table  2 

-5 

(pg.92).In  every  case  F(x)  +  F(-x)  =  1+e  where  |e|  <  5*10  .  The 


n  ■  ■  :  ...  '• " 
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apparently  bad  solutions  for  a  =  0.9  and  0.95  do  not  indicate  a  failing 
of  the  procedure,  but  merely  reflect  the  fact  that  as  a  -> 1  ,  the  arguments 
(x+l-a^)/a^  and  (x-l4a0)/a2  tend  to  x  so  that  the  functional  equations 
(4.1)  become  a  tautology.  By  halving  the  tabular  interval  so  that  h  =  0.05, 
the  resulting  solutions  appear  to  be  accurate  to  six  or  seven  figures  (cf. 
pages  101  and  102). 

Another  check  on  the  procedure  is  the  comparison  of  the  numerical 
solution  by  a  system  of  linear  equations  with  a  Monte-Carlo  solution  given 
by  Bush  and  Mosteller  in  [1]  p.133*  These  authors  calculated  on  approximate 
limiting  distribution  by  a  1000  -  trial  stat-rat  for  the  parameter  values: 

*1  =  *2  =  0,50  ’ 

=  0 . 60  ,  oig  =  0 . 90  ; 

'hl  =  0.75  ,  a2  =  0.10  . 

The  results  of  their  calculation  are  displayed  graphically  in  figure  6. 

Since  the  numerical  solution  by  the  method  of  linear  equations  tabulates  the 
transformed  distribution,  an  inverse  transformation  was  made  to  map  the 
limiting  distribution  onto  the  domain  [  .10,  .75],  The  results  of  this  solution 
are  displayed  graphically  in  figure  7 . 
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Figure  6  ; 


Figure  7 


F(p) 


0  0.2  0.4  0.6  0.8  1.0 


The  approximate  asymptotic  cumulative  distribution  of  p  obtained 
from  a  1000-trial  stat-rat  for  Experimenter-Controlled  events 
Q^p  =  0.3  +  0.6  p,  Q2p  =  0.01  +  0.9  p,  ^  =  tt2  =  0.5  .  The 

vertical  dotted  lines  show  the  trapping  limits  X2  =  0,10  and 

\  =  0.75.  (Reproduced  from  [1]  p.  135). 


0  0.2  0.4  0.6  0.8  1.0 


The  approximate  limiting  distribution  of  p  obtained  from  a 
system  of  39  linear  equations  in  F(p)  for 

^=*2=0.5  ,  a1=0.60  ,  a2  =  0.90,  0.75  ,  ^2=0.10  . 


•  -  ' 

fr 
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A  second  set  of  tables  was  calculated  for  values  of  the  parameters 
which  yield  known  absolutely  continuous  distributions  (Zidek  [9])*  These 
tables  include  the  parameter  values  ,  a.j=ap=i^n  ,  (n=2,  3,4,  ,  . . ,  10) . 

A  number  of  spot  checks  were  made  by  calculating  numerical  values  of  the 
distribution  from  the  known  absolutely  continuous  functions  and  comparing 
these  results  with  values  taken  from  the  numerical  tables.  Table  1 
displays  some  of  these  checks.  Since  the  functions  are  centro-symmetric, 
only  values  of  F(x)  for  x  e(-l,  0)  need  be  included.  Note  that  the 
maximum  error  is  <  +  2x10  "r  . 

A  final  set  of  tables  was  calculated  for  some  representative 
values  of  the  parameters  to  display  the  behaviour  of  the  distributions 
for  a^+a2  —  ^  ' 

Although  it  has  not  been  proved  that  the  numerical  procedure 
converges  to  the  true  solution,  the  above  tests  indicate  that  a  reasonable 
amount  of  confidence  may  be  placed  on  the  numerical  solution. 
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n 

X 

True  Value  of 

F(x) 

Approximate  Value  of 

F(x) 

2 

-.900 

.005  151 

.005  151 

-.750 

.032  197 

.032  197 

-  .600 

.082  430 

.082  430 

-.450 

.155  835 

.155  850 

0 

0 

KN 

1 

.252  426 

.252  431 

-.150 

.371  960 

•371  957 

3 

-.900 

.000  593 

.000  593 

-.800 

.004  746 

.001  746 

O 

O 

8 

.127  834 

.127  828 

-.200 

.288  896 

.288  868 

-.100 

.390  917 

.390  880 

4 

0 

0 

ON 

1 

.000  072 

.000  069 

-.700 

.005  817 

.005  817 

-.500 

.044  871 

.044  740 

-.4oo 

.090  959 

.091  001 

-.200 

.256  953 

.256  741 

6 

1 

VO 

0 

0 

.000  oil 

.000  000 

-.800 

.000  073 

.000  068 

-.700 

.000  829 

.000  824 

Table  1. 


Comparison  of  Values  of  the  distribution  function 
calculated  from  the  absolutely  continuous  functions  and 
those  from  the  numerical  tables  for  jt. 
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ALPHA  1 =  .60000000  P  1 1  = 
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ALPHA 1 =  .9^000000  PI1= 
X 

-1.000000 
-.950000 
-.900000 
-.850000 
-• 800000 
-.750000 
-.700000 
-.650000 
-.600000 
-.550000 
-.500000 
-.450000 
-.400000 
-.350000 
-.300000 
-.250000 
-.200000 
-.150000 
-.100000 
-.050000 
0.000000 
.050000 
.100000 
.150000 
.200000 
.250000 
.300000 
. 3^0000 
.400000 
.450000 
.500000 
.550000 
.600000 
.650000 
.700000 
.750000 
.800000 
.850000 
.900000 
.950000 
1.000000 


500  ALPHA2=  .95000000 

F 

0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
0.000000 
. 000019 
. 000125 
.000547 
. 001883 
.005436 
.013594 
.030062 
.059641 
.107338 
.176811 
.268617 
.379006 
.500000 
.620993 
.731382 
.823187 
.892660 
.940357 
. 969937 
.986404 
. 994563 
.998116 
. 999452 
. 999874 
. 999980 
. 999999 
1.000000 
1.000000 
1.000000 
1.000000 
1 . 000000 
1.000000 
1 . 000000 
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Table  5  below  displays  the  numerical  solutions  of  the  absolutely 

\  1  /n 

continuous  distribution  functions  corresponding  to  3t^=Ttp='?  »  a^=a^=  (■!■)  ' 

(n=2, • • • > 10 ) .  The  following  tables  have  been  rounded  to  four  significant 
digits . 
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ALPHA 1 =  .70710678  P 1 1  = 
X 

-1.000 
-.950 
-.900 
-.850 
-.800 
-.750 
-.700 
-.650 
-.600 
-.550 
-.500 
-.450 
-.400 
-.350 
-.300 
-.250 
-.200 
-.150 
-.100 
-.050 
0.000 
.050 
.100 
.150 
.200 
.250 
.300 
.350 
•  4  00 
•  4^0 
.500 
.550 
.600 
.650 
.700 
.750 
.800 
.850 
.900 
.950 
1.000 


.500  ALPHA2=  .70710678 

F 

0.0000 
.0013 
.0052 
.0116 
.  0206 
.0322 
.0464 
.0631 
.0824 
.1043 
.1288 
.1559 
.1855 
.2177 
.  2524 
.2898 
.3297 
.3720 
.4146 
.4573 
.  5000 
.5427 
.  5854 
.6280 
.6703 
.7102 
.7476 
.  7823 
.8145 
.  8441 
.8712 
.  8957 
.9176 
.9369 
.9536 
.9678 
.9794 
.9884 
.  9948 
.  9987 

1.0000 


TABLE  3.  NUMERICAL  SOLUTION  OF  THE  ABSOLUTELY  CONTINUOUS  LIM 

DISTRIBUTIONS  FOR  N=2-10. 
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ALPHA  1 =  .79370053  PI1  = 
X 

-1.000 

-.950 

-.900 

-.850 

-.800 

-.750 

-.700 

-.650 

-.600 

-.550 

-.500 

-.450 

-.400 

-.350 

-.300 

-.250 

-.200 

-.150 

-.100 

-.050 

0.000 

.050 

.100 

.150 

.200 

.250 

.300 

.350 

.400 

.450 

.500 

.550 

.600 

.650 

.700 

.750 

.800 

.850 

.900 

.950 

1.000 


.500  ALPHA2=  .79370053 

F 

0.0000 
.0001 
.  0006 
.0020 
.  0047 
.0093 
.0160 
.0254 
.0380 
.0541 
.0742 
.  0987 
.1278 
.1616 
.1999 
.2425 
.2889 
.3385 
.  3909 
.4450 
.5000 
.  5550 
.6091 
.6615 
.7111 
.7575 
.  8001 
.  8384 
.  8722 
.9013 
.9258 
.  9459 
.9620 
.  9746 
.  9840 
.9907 
.9953 
.  9980 
.  9994 
.9999 

1.0000 
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c  " ^  A  H  P  ..J  A 

r  ~ "  #  =  r  i  o 

.  r  c  n  -  ' 

x 

r-  r\  r  r* 

<* 

r  r>  ■> 

• 

0  -  p  .  - 

a  o  n  r  . 

r  r  _ 

• 

0^00  . 

TJ  p  _ 

P  A  ~  r  . 

Q  .- 

■  o  r.  r  „ 

"  3r  _ 

or  r . 

+-  0  C  . 

-3a,-- 

r.  o  r 

r  A  3  ^  . 

■a - .  - 

0  A  v'  . 

00 

TQOP  ^ 

a  A  .  - 

C,  T  c  r  . 

' r '  r . 

C  3  0  .  _ 

o  o  o  r 

0  0.- 

. 

0  n  0  •  - 

O  3  q  C  . 

00  0 

?  Q  <  '  . 

' 3  r .  - 

P  P  C  r  # 

r  .  _ 

4 

■'  o  r  ~>  . 

0  0.0 

0  3  f  3  # 

0  '  '  . 

r  O  0  r>  . 

r\  r  r 

-  r  a  . 

0  3  r . 

r  r  , 

0  r 

3  —  3  r 

• 

r  1  r 

r  r  o  o  . 

-  r  r- 

A  Q  °  p  . 

3  " 

00r'  . 

00  A. 

c  I  '  -■ . 

0  3  A  . 

°  =  S  "  . 

r  ,a  ^  o  . 

0  3  3 

'  A  °  . 

>  . 

a ''To . 

a. 

0  a  «  ? . 

•  />  r- 

• 

r  p .  o  r 

rj 

r  p  o  r>  # 

0.  . 

0 

3  Q  . 

Aoor  . 

^  p 

o  o  on 

3P  . 

0 00'  .  r 

00  '  .  r 

r,  r  *HC  J,» 


r  1  j  p  a  T 


106- 


ALPHA 1 =  • 8408  9642  P  1 1  = 
X 

-1.000 
-.950 
-.900 
-.850 
-.800 
-.750 
-.700 
-.650 
-.600 
-.550 
-.500 
-.450 
-.400 
-.350 
-.300 
-.250 
-.2  00 
-.150 
-.100 
-.050 
0.000 
.050 
.100 
.150 
.200 
.250 
.300 
.350 
.400 
.450 
.500 
.550 
.600 
.650 
.700 
.750 
.800 
.850 
.900 
.950 
1.000 


.500  ALPHA2=  .84089642 

F 

0.0000 

0.0000 
.0001 
.0004 
.0011 
.0028 
.0058 
.0108 
.0184 
.0294 
.0447 
.0650 
.0910 
.1231 
.1615 
.2062 
.2567 
.3125 
.  3725 
.4355 
.5000 
.5645 
.6275 
.6875 
.  7433 
.  7938 
.  8385 
.8769 
.  9090 
.9350 
.  9553 
.  9706 
.9816 
.  9892 
.9942 
.9972 
.  9989 
.9996 
.  9999 

1.0000 

1 . 0000 
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ALPHA 1 =  .87055056  P 1 1  = 
X 

-1.000 
-.950 
-.900 
-.850 
-.800 
-.750 
-.700 
-.650 
-.600 
-.550 
-.500 
-.450 
-.400 
-.350 
-.300 
-.250 
-.200 
-.150 
-.100 
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.850 
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1.000 


.500  ALPHA2=  .87055056 

F 
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0.0000 
0.0000 
.0001 
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.  0047 
.  0092 
.0164 
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ALPHA  1 =  .89089872  P  1 1 
X 

-1.000 

-.950 

-.900 

-.850 

-.800 

-.750 

-.700 

-.650 

-.600 

-.550 

-.500 

-.450 

-.400 

-.350 

-.300 

-.250 

-.200 

-.150 

-.100 

-.050 

0.000 

.050 

.100 

.150 

.200 

.250 

.300 

.350 

.400 

.450 

.500 

.550 

.600 

.650 

.700 

.750 

.800 

.850 

.900 

.950 

1.000 


.500  ALPHA2=  .890898 

F 

0.0000 
0.0000 
0.0000 
0.0000 
.0001 
.  0003 
.0008 
.0021 
.0046 
.0093 
.0172 
.0297 
.0483 
.0744 
.1093 
.1538 
.2081 
.2716 
.3429 
.4200 
.5000 
.  5800 
.6571 
.7284 
.7919 
.  8462 
.  8907 
.9256 
.9517 
.9703 
.9828 
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.  9999 
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1 .0000 
1 . 0000 
1.0000 
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ALPHA 1 =  .90^72366  P  T 1  = 
X 

-1,000 

-.950 

-.900 

-.850 

-.800 

-.750 
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-.650 

-.600 

-.550 

-.500 

-.450 

-.400 

-.350 

-.300 

-.250 

-.200 

-.150 

-.100 

-.050 

0.000 

.050 

.100 

.150 

.200 

.250 

.300 

.350 

.400 

.450 

.500 

.550 

.600 

.650 

.700 

.750 

.800 

.850 

.900 

.950 

1.000 


.^00  ALPHA2=  .90572366 

F 

0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
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.0009 
.0024 
.0053 
.0109 
.0204 
.0357 
.  0587 
.0910 
.1342 
.1889 
.2548 
.3304 
.4133 
.  5000 
.5867 
.6696 
.  7452 
.8111 
.  8658 
.  9090 
.  9413 
.9643 
.9796 
.  9891 
.  9947 
.9976 
•  99P 1 
.  9997 
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1 . 0000 
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1 . 0000 
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1 . 0000 
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The  following  table  displays  the  behaviour  of  the  transformed 
limiting  distribution  functions  for  selected  values  of  the  parameters. 

Figure  8  page  113  contains  a  graph  of  the  absolutely  continuous 
distributions  for  a  =  (-4)^^  (see  also  the  corresponding 
results  of  Zidek  [9])- 
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CHAPTER  V 

NUMERICAL  TABULATION  OF  THE  TRANSFORMED  LIMITING 
DISTRIBUTIONS  OF  THE  SUBJECT  CONTROLLED  MODEL 

§  5. 1  Introduction 

Because  of  the  difficulties  involved  in  the  Subject-Controlled 
Model,  numerical  solutions  may  be  the  only  possible  method  of  finding  the 
forms  of  the  limiting  distributions.  Exact  solutions  of  the  known  purely 
singular  cases  of  <  1  cannot  be  constructed  in  general,  and  no 

absolutely  continuous  distributions  have  been  found. 

However,  the  functional  integral  equations  satisfied  by  the 
limiting  distributions  may  be  solved  numerically  to  obtain  approximations 
of  the  limiting  distribution  functions  for  all  values  of  the  parameters. 
The  accuracy  of  the  method  cannot  be  theoretically  determined,  but  in 
the  centro- symmetric  cases  the  numerical  solutions  appear  to  be 
sufficiently  accurate  for  practical  purposes.  Comparing  the  numerical 
solutions  with  Monte-Carlo  approximations  also  indicates  that  the 
numerical  tabulation  is  reliable  and  it  can  be  used  to  obtain  greater 
precision  than  is  feasible  by  the  Monte-Carlo  method, 

§5.2  Numerical  Procedure 

Consider  the  problem  of  calculating  the  values  of  the 
distribution  functions  from  the  functional  integral  equations  on  a 
tabular  set  of  points  (x^)  lying  in  the  domain  (-1,1).  These  points 
are  chosen  such  that  x^  =  ih-1  (i  =  0,  1,  2,  ...,  2N)  where  h  =  1/fo. 


- 
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On  this  tabular  set,  the  functional  integral  equations  may  be  written 

(5.1)  -F(x. )+[*.+( x  +l)(Ap-A1)/(2a .)]F[(x  +l-a  )/a  ] 

(x^+l-a^  )/a^ 

+[l-A2+(l-x1)(>l2->1)/(2a2)]F[(x1-l+a2)/a2]-|^2->v1)  J  F(u)du  =  0  , 

(x.-l+a2)/a2 

<*2  >  \), 

and 

(5.2)  -F(x  )+[A  +(x.-l)(A  -A2)/(2a1)]F[(xi-l+a1)/a1] 

(xi+l-a2VCL2 

+[l-X2-(xi-l)(A1-X2)/(2a2)]F[(x1+l-a2)/a2]4(>1->21  J  F(u)du  =  0  , 

(xi-l+a1)/a1 

(\  >\,  \/l,  0  )  j 

where  (i  *  1,  2,  3,  2N-l). 


Imposing  the  lower  boundary  condition  of  F(x)  =0,  x  <  -1  , 
the  above  equations  become 

(x.+l-c^  )/ax 

(5.3)  -F(xi)+[X1+(x1+l)(7v2-X1)/(2a1)]F[(x1+l-a1)/a1]4(A2-X1)  j  F(u)du  =  0  , 

-1 

(i  =  1  j  2,  5,  ...  m)  j 

where  mh-1  <  l-2a^  ,  (m+l)h-!  l-2a^  , 

and 

(x.+l-a2)/a2 

(?A)  -F(x1)+[l-A2-(x1+l)(X1->v2)/(2a2)]F[(x1+l-a2)/a2]+i(X1-?l2)  j  F(u)du  =  0, 

-1 


(^  —  2,  3;  •  •  •  > 


m), 


.  ’(l*>  -  .  A  -  l4„  ) 
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where  mh-1  <  l-2a^  ,  (tn+l)h-l  >  l-2a^ 

Imposing  the  upper  boundary  condition  gives 

(5.5)  -F(x1)+A2+[l-X2-(l-x.)(X2-A1)/(2a2)]F[(x1-l+aa)/a2] 

1 

-j(VV  J  F(u)du  =  0  ,  (1  =  n,  n+1,  2N-1)  , 

(Xi'1+a2>/a2 

where  (n-l)h-l<  2a^-l  ,  nh-1  >  2a^-l  , 

and 

(5.6)  -F(xi)+l-^1+[X1+(xi-l)(X1-X2)/(2a1')]F[(xi-l+a1)/a1] 

1 

+  2  J  F(u)du  =0  ,  (i=n,n+l,n+2, . . . ,2N-l)  ; 

(xi-l+a1)/a1 

where  (n-l)h-l  <  2a^-l  ,  nh-1  >  2a^-l  . 

As  in  Chapter  IV,  the  arguments  (x^+l-a^ )/a-^  ,  (x^-I-hx^)/0^  > 
(x^+l-a^ )/a'2  >  (xjL“-*-+ai)/ai  will  not,  in  general,  lie  on  the  tabular 
points  (xi).  Thus  equations  (5.1),  (5.3),  (5.5)  and  (5.2),  (5.4),  (5.6) 
define  systems  of  2N-1  equations  in  6N-m-n~5  unknown  values  of  F 
as  well  as  2N-1  undetermined  integrals.  The  values  of  the  distribution 
function  at  off-tabular  points  may  be  expressed  in  terms  of  its  values  on 
adjacent  tabular  points  by  the  Lagrange  Interpolation  formula  described 
in  Chapter  IV, 

The  integrals  in  each  of  the  equations  (5.l)-(5.6)  can  be 
approximated  numerically  by  weighted  sums  of  F(x^)  ,  F  [(x^+l-a1  )//a'  ] 


>nr-  ■  .3-  ■  :  ...  t-»aaif5 
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and  F [ (x^-l+a* )/o>*]  ,  where  a'  and  a*  are  to  be  interpreted  as 
shown  on  page  71  .  Wherever  possible,  a  Simpson's  Rule  approximation 
was  used,  but  many  situations  arose  where  the  simpler  Trapezoidal  sum 
was  required.  The  following  notation  enables  us  to  discuss  a  number  of 
cases  together. 

Let  the  largest  argument  be  (x+l-a'  )/&'  and  the  smallest 
be  (x-l+a*)^*.  For  x  =  x^  ,  let  (x^+l-a*  )/a»'  be  denoted  by 
and  (x^-l+a* J/a*  be  denoted  by  , 

Suppose  x.  ,  <  A.  <  x,  and  x.  ,  <  c.  <  x. 

k~  1  i  —  k  j"l  i  —  j 


Figure  9  :  Arrangement  of  the  arguments  A^  =  (x^+l-a'  )/cl'  and 

=  (x^-l+a*)^*  when  x  =  x^ 

The  integral  to  be  approximated  is  of  the  form 

A. 

x 

(5.7)  f  F(u)  du  . 

ci 

Each  of  the  following  cases  requires  special  consideration: 


■■■  ■ 
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Note! 


Ai  /  xk 
i)  k-j  <  2 

F(u)du  ~  l(xj-Ci)[F(Ci)+F(xi)]+i(Ai-xi)[F(xi)+F(Ai)] 


A. 

n  1 


ii)  k-j  =  2 
A. 


J 

C. 

i 


/  F(u)du  -v  i(x  -C.  )[F('C.  )+F(x  )]+Jh[F(x  )+F(x  .  )J 

i  “  J  i  1  J  J  J+i 

+l(Ai-xk_1)[F(xk_1)+F(Ai)]  . 


iii )  k-j  =  2 v- 1  ,  r  =  2,  3,  4,  . . . 


F(u)du  ~  i(x  -C.)[F(C.)+F(x  )] 
~  J  L  J 


(k-J4)/2  (h-t?)/6 

+(V5)'F(xj)^  F(xj+2J!.l)  +  2  )  F(xj+2i)+F(xk-l 

i=l  i=l 


+  i(Ai-xk_1)[F(xk_1)+F(AjL)3  . 


iv)  k-j  =  2r  ,  r  =  2,  3,  4,  . . . 


+(h/3r  F(x 

L 

+ih[F(xk_1 


C  )[F(c.)+F(xj)]+ih[F(xj)+F(xj+1)] 

(k-J-2)/2  (k-J-4)/2 

j+1^4  ),  F^Xj+l+2i)  +  2  F(xj+2+2i 

i=l  i=l 

)+F(Ai)]. 


)+F(xk 


when  C  ^ 
when  A^ 
and  when 


<  -1  > 
>  1  > 


j  =  1  and  equation  (5.3)  or  (5.4)  is  used, 
k  =  2N  and  equation  (5.5)  or  (5.6)  is  used, 
the  first  Trapezoidal  Sum  is  vacuuous. 
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b)  \  = 

i)  k-j  <  2 
i 

F(u)du  ro  i(xi-Ci)[F(Ci)+F(xi)]  +  i(h)[F(xi)+F(xi+1)] 


ii)  k-j=2 
A 


F(u)du  =  ^(x  -C  ) [F(C  )+F(x  ) ]+(h/3) [F(x  )+4F(x  )+F(x  ) ] 


J  i  i  J 


j  +  1  j+2‘ 


iii)  k-j=2r-l  r  =  2,3,1)-,... 
A 


iv) 


f  F(u)du  -v  i(xj-C1)[F(C1)+F(xj)]+|h[F(xj)+F(xj+1)] 

1  |(k- j-1) 

(h/3)[F(x.+1),  4  /  F(x.+2)+  2  F(x.+  1+2<) 

L  4=1  4=1 


+  F(xk) 


k- j=2r  r  =  2,3, . 

A 


f  1  F(u)du  -v  Kx.-C^tFtC.KFtx.)] 


c. 

1 


l(k-i-u) 

+(h/3)  F(x.)+  1*  )  F(xJ+1+2  )+  2  ^  F(x.+2+2<)+F(xk) 

L  4=0  4=0 


When  Ch  <  -1,  j=l  and  (5.3)  or  (5.^)  is  used, 

when  >  1 ,  k=2N  and  (5-5)  or  (5.6)  is  used^ 

when  =  x^ ,  the  first  trapezoidal  sum  becomes  vacuuous. 

In  all  cases,  when  or  A^,  does  not  lie  on  a  tabular  point, 

F(Ch)  and/or  F(A^)  is  expressed  in  terms  of  the  Lagrange  Interpolation 


polynomial . 
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The  ALGOL  form  of  the  program  for  setting  up  the  system  of  linear 
equations  follows. 


procedure 

Subject-Controlled  Distribution  Functions; 

array 

A  [  1 : 50 ,1:51],  x[0:51],  y[l:50],  C[l:4]; 

real 

h,  alfal,  alfa2,  laml,  lam2,  argl,  arg2,  coefl,  coef2,  coef3, 

swtch,  all,  al2,  B,  sgn,  sgnl,  hf,  ht,  d,  e  ,  comax,  comin,  arg, 

cof,  hp,  hpstor; 

integer 

i,  j,  k,  max,  max2,  maxi,  lead,  init,  jstor,  kick,  no; 

Boolean 

kode,  jump,  mop,  leap; 

begin  comment 

This  procedure  tabulates  the  values  of  the  limiting  distribu¬ 
tion  functions  by  calculating  the  numerical  solutions  of  the 

functional  integral  equations  satisfied  by  the  limiting  dis¬ 
tributions  of  transformed  response  probabilities; 

read  (h,  alfal,  alfa2,  laml,  lam2); 

max:=  2.0/h;  max  2:=  max+1; 

for  i:=0  step  1  until  max  2  do  xfi]:=i  X  h  -  1.0; 

if  laml  >  lam2  then  go  to  type  2  else  go  to  type  1; 

type  1: 

kode:=true;  argl;=1.0-alfal;  arg2: =alfa2-one; 

coefl:=lam2-laml;  coef2: =coef l/(2.C X  alfal); 

coef 3 :=coef 1/(2 .0  x alfa2) ; 

swtch: =1 .0-2.0 Xalfa2;  all;=alfal,  al2:=alfa2;  B:=-lam2; 

sgn:=  -1.00  ; 

go  to  calculate: 

type  2: 

kode:=false:  arel:=1.0-alfa2;  are2:=alfal-1.0:  coef 1: =laml- lam2 

coef2:=coef l/(2.C  Xalfa2);  coef3:=coefl/(2.0Xalfal); 

swtch;=l .0-2.0Xalfal;  all:=alfa2;  al2:=alfal;  B:=laml-1.0; 

sgn: =1.0  ; 

-  ■ 


A-  f 

■ 


' 


'  '  ' 
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comment 


calculate: 


begin  comment 


comment 


The  above  program  segment  reads  in  the  parameters  and  initial¬ 
izes  for  the  required  functional  integral  equation; 
maxi: =max-l;  hf : =2.0XhXcoef lXsgn/3. 000; 
ht:=hX-oef lXsgn/3.000;  sgnl:=sgn;  leap; =true; 
d: =hxsgnXcoef lx. 25000 ;  e: =sgnxcoef lx, 250000; 
for  i=l  step  1  until  max  1  do 

This  loop  sets  up  the  i' ^  row  of  the  matrix  A[i,j]; 
for  j:=l  step  1  until  max  do  A[i,j]:=0.0; 

A[i, 1]:=-1.0 

if  kode  then  begin  comax:=(x[i]+1.0)xcoef2+laml; 

comin:=( 1 ,0-x[i])xcoef3+1.0-lam2; 

end 

else  begin  comax:=l,0-lam2-(x[i]+1.0)xcoef2; 

comin: =(x[ i ] - 1 . 0) xcoef 3+laml ; 

end 

if  leap  then  go  to  pace  else  go  to  run; 

pace:  arg: =(x[i]+arg l)/all;  coef:=comax; 

if  arg  <  1.0  then  go  to  ii.ci.ch  else  go  to  trot; 

march:  for  j:=l,  j+1,  while  x[j]  <  arg  do; 

if  x[j]  =  arg  then  go  to  equal  else  go  to  unequal; 

equal:  hp:=0,0;  A[i, j]:=A[i, j]+cof ;  go  to  test; 

unequal:  Lf  jump  then  hp:=arg-x[j]  else  hp: =x[ j ] -arg; 

_if  j  <  3  then  lead;=0  else 

if  1  >  maxi  then  lead:=max-4  else  lead:=j-4; 

Interp  (lead,  max,  arg,x,C); 

Interp  is  the  Lagrange  Interpolation  routine  described  in 


Chapter  IV; 


:s  ^ 
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trot: 

run: 

test : 

skip: 


tabpt: 


nontabpt : 


set: 


for  k  ;  =  1  step  1  until  4  de  A[  i  ,k+lead  ] :  =A[  i  ,k+lead  ]+c(k  )x(coefH-hj  x,e  ) ; 
go  to  test; 

leap:=false;  sgnl:=-sgn; 

hpstor:=0.0;  A[i,max]:=A[i,max]4B ;  jstor:=max; 

go  to  skip' 

if  jump  then  begin  jstor:=j;  hpstor:=hp;  go  to  skip;  end 
else  go  to  set; 

if  x[i]  >  swtch  then  begin  cof:=comin;  arg: =(x[ i ]+arg2)/al2; 
jump: =false;  go  to  unequal;  end 
else  begin 

A[ i, jstor-1 ] : =A[ i, jstor-1 ]+hpstorxe;  1 :  =  j  stor ; 
if  hpstor=0  then  go  to  tabpt  else  go  to  nontabpt; 
end 

if  l/2.0-entier  (l/2.0)=0  then  begin 

kick: r 1;  no: = jstor-1;  go  to  integral  end 
else  begin  A[ i , 1 ] : =A[ i, 1 ]+d;  kick:  =2  ;  no:=jstor-2  go  to  stintegral 
end 

if  l/2„0-entier  (l/2.0)=0  then  begin 

A[i,l]:=A[i,l]+d;  kick;=2;  no:=jstor-2;  go  to  stintegral* 
end 

else  begin  kick:=l;  no:=jstor-2;  g°  to  integral;  end 
A[i, jstor-1]: =A[ i , jstor-1 ]+Hps tor xe;  A[i,j]:=A[i,j ]+hpxe; 
l:=jstor-j; 

if  hpstor=0  then  go  to  tabpt  2  else 
begin  if  1  <  2  then  go  to  main  else 
if  1=2  then  go  to  place  else 
begin  if  1/2.0-entier  (l/2.0)=0  then 


: 

.  ‘  ‘  J-V’  .;  ;u^L 

-■■  ......  .  - 

— ■'  ■*»—  -•■MW 


-134- 


begin  kick:=j+2;  A[i, j ] : =A[i, j ]+d; 

A[ i , j+1 ] : =A[ i , j+1 ]+d;  no:=jstor-2; 
go  to  stintegral; 

end 

else  begin  kick:=j+l;  no:=jstor-2;  go  to  integral*  end 

end 

end 

tabpt  2:  if  1  <  2  then  go  to  place  else 

begin  if  1=2  then  begin  kick:=j+l;  no:=jstor-l; 

go  to  stintegral;  end 
else  if  l/2»0-entier( l/2.0)=0  then 
begin  kick:=j+l;  no:=jstor-l; 
go  to  stintegral; 

end 

else  begin  kick:=j+2;  A[i, j ]:=A[i, j ]+d; 

A[ i, j+1 ] : =A[ i, j+1 ]+d;  no;=jstor-l; 
go  to  stintegral; 

end 

place:  A[i,j]:=A[i,j  ]+d;  A[ i,  j+1  ] :  =A[  i,  j+1  ]+sgnlxhx:coef  lx.  25000; 

go  to  main; 

stintegral:  A[ i, kick -1 ]=A[i, kick- 1 ]+sgnxhxcoef 1/6.0000; 
integral:  mop: =true * 

for  k:=kick  step  1  until  no  do 

begin  if  mop  then  begin  A[i,k]:=A[i,k]+hf ;  mop: =false;  end 
else  begin  A[i,k]:=A[i,k]+ht;  mop: =true;  end 

end 

A[  i  ,no  ] :  =A[  i  ,no  ]+sgn  Ixb  xcoef  1/6.000; 


" 

- 

•; 

. 

. 
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main:  end  This  ends  the  loop  for  setting  up  the  matrix  equation; 

GAUSEL  (maxi,  A,  y);  comment  This  statement  calls  a  Gaussian 

elimination  with  interchanges  routine  to 
solve  the  matrix  equation; 

Punch  (laml,  alfal);  Punch  (lam2,  alfa2);  Punch  (-1,0); 

for  i=l  step  1  until  maxi  do  Punch  (x[i],y[i]); 

Punch  ( 1 , 1 ) ; 

end  Subject-Controlled  Distribution  Functions; 

§3.5  Discussion  of  the  Numerical  Results. 

The  procedure  described  above  was  programmed  in  the  Fortran  II 
language  and  several  examples  were  run  on  the  University  of  Alberta's 
IBM  1620  computer.  Hardware  floating  point  with  a  seven  decimal  digit 
maitissa  and  two  decimal  digit  exponent  was  used.  The  parameter  values 
chosen  were  the  ones  used  by  Bush  and  Mosteller  in  calculating  approximate 
distributions  by  Monte-Carlo  Methods, 

Useful  error  bounds  cannot  be  obtained  in  the  Subject  case  since 
practically  nothing  is  known  of  the  properties  of  the  limiting  distributions 
We  cannot  even  compare  the  numerical  values  with  known  solutions.  However, 
there  is  one  check  that  indicates  the  numerical  procedure  is  reliable.  This 
is  given  by  the  centro- symmetric  distributions  of  ^i+^2=^  ([l)pl07) 

Two  of  the  numerical  examples  satisfy  this  condition  and  with  a  tabular 
interval  h-.05,  the  error  in  centro -symmetry  is  <  2x10”^.  (see  Table  5) 

Comparison  of  the  results  of  the  Monte-Carlo  approximations  with 
the  numerical  solutions  of  the  functional  integral  equations  indicates  that 
the  Monte-Carlo  results  tend  to  be  much  smoother.  Since  the  Monte-Carlo 
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approximations  were  determined  by  at  most  four  thousand  trials,  the  maximum 
accuracy  possible  is  one  to  two  decimal  digits.  On  the  other  hand,  the 
numerical  solution  of  the  functional  integral  equations  seems  to  yield  an 
accuracy  of  four  to  five  digits. 
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LAMBDA1 =  0.00000 
LAMBDA2=  1.00000 


X 

-1.0000 

-.9500 

-.9000 

-.8500 

-.P000 

-.7500 

-.7000 

-.6500 

-.6000 

-.5500 

-.5000 

-.4500 

-.4000 

-.3500 

-.3000 

-.2500 

-.2000 

-.1500 

-.1000 

-.0500 

0.0000 

.0500 

.1000 

.1500 

.2000 

.2500 

.3000 

.3500 

.4000 

.4500 

.5000 

.5500 

.6000 

.6500 

.7000 

.7500 

.8000 

.8500 

.9000 

.0500 

1.0000 


ALPHA1 =  .500000 
ALPHA2=  .500000 


F 

0.0000 

.0001 

.0013 
.0044 
.0164 
.0249 
.0404 
.0841 
.1102 
.1431 
.  1500 
.1581 
.2040 
.2478 
.3364 
.3749 
.4062 
.4641 
.4850 
.4987 
.5000 
.5013 
.5150 
.5359 
.5938 
.6251 
.6636 
.7522 
.7960 
.8419 
.8500 
.8569 
.8898 
.9159 

.9596 

.9751 

.9836 

.9956 

.9987 

.9909 

1.0000 


TABLE  5.  NUMERICAL  APPROXIMATION  OE  SOME  LIMITING  0 
THE  SUBJECT-CONTROLLED  MODEL  EOR  SELECTED  VALUES  OF 


STRIBUTIONS  OE 
THE  PARAMETERS 
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L  AMRn A 1 =  , ft  0000 

LAMBDA2=  .20000 


X 

-1.0000 
-.9500 
-.9000 
-.8500 
-.8000 
-.7500 
-.7000 
-.6500 
-.6000 
-.  5500 
-.5000 
-.4500 
-.4000 
-.3500 
-.3000 
-.2500 
-.2000 
-.1500 
-.1000 
-.0500 
0.0000 
.0500 
.1000 
•  1500 
.2000 
.2500 
.3000 
.3500 
.4000 
.4500 
.5000 
.5500 
.6000 

.6500 

.7000 

.7500 

.8000 

.8500 

.9000 

.9500 

1.0000 


ALPHA  1  =  .500000 
ALPHA2=  .500000 


F 

0.0000 
.1484 
.1877 
.2219 
.2395 
•  2  644 
.2866 
.2974 
.3120 
.3232 
.3496 
.3750 
.3845 
.3952 
.4022 
.4147 
.4275 
.4353 
.4482 
.4607 
.5000 
.  5393 
.5518 
.5647 
.5725 
.  5852 
.5978 
.6048 
.6155 
.6250 
.6504 
.6768 
.6880 
.7026 
.7133 
.7356 
.7605 
.7781 
.8123 
.8516 
1.0000 
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LAMBDA  1  = 

.75000 

ALPHA1 = 

.600000 

L AMBDA2  = 

.20000 
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.700000 
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CHAPTER  VI 

CONCLUSION 

The  functional  equations  satisfied  by  the  limiting  distribution 
functions  of  the  Bush  and  Mosteller  stochastic  learning  models  have  yielded 
much  valuable  information  about  the  limiting  distribution  of  response 
probabilities.  Exact  solutions  for  the  distribution  functions  can  be 
derived  from  these  equations  in  many  cases  and  numerical  solutions  of  the 
functional  equations  appear  to  be  a  fast  and  accurate  method  of  tabulating 
the  limiting  distributions  for  all  values  of  the  parameters. 

The  distribution  functional  equation  derived  from  the  Experimenter- 
Controlled  Model  has  been  particularly  valuable.  For  all  cg^-kx^  <  1  ,  it 
can  be  used  to  construct  the  limiting  distribution  function  by  a  very  simple 
mechanical  algorithm.  For  >  1  ,  the  functional  equation  yields  a 

construction  for  absolutely  continuous  distribution  functions  where  they 
exist.  By  pursuing  this  approach  further,  more  precise  criteria  for  absolute 
continuity  might  be  derived.  However,  even  where  analytic  methods  fail, 
numerical  solutions  appear  to  provide  sufficient  information  to  be  of  value 
for  applied  problems. 

The  distribution  functional  integral  equations  of  the  Subject- 
Controlled  Model  have  not  yet  proved  to  be  of  particular  value  For  a^+a^  <  1* 
the  functional  equations  merely  demonstrate  the  form  of  the  limiting  distributions. 
However,  future  work  may  yield  a  method  of  constructing  exact  solutions  from 
these  apparently  recalcitrant  equations.  As  for  the  distribution  functions 
corresponding  to  values  of  a^-KX^  >  1  ,  almost  no  useful  information  has 
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been  derived.  Despite  the  poverty  of  analytical  results,  the  functional 
integral  equations  do  appear  to  be  amenable  to  numerical  methods.  Even 
when  (where  uniqueness  of  the  solutions  has  not  been  established) 

the  numerical  results  seem  to  be  stable.  This  tends  to  imply  that  the 
solutions  are  unique. 

Until  more  is  learned  about  the  structure  of  these  types  of 
functional  equations,  we  may  have  to  depend  upon  numerical  results.  It 
may  even  be  hoped  that  numerical  work  can  serve  as  a  guide  to  theoreticians 


in  this  field. 


-  143  - 


BIBLIOGRAPHY 


[1]  Bush,  R.R.  and  Mosteller,  F:  Stochastic  Models  for  Learning.  Wiley 

(1955). 

[2]  McGregor,  J.R.  and  Hui,  Y.Y;  ’’Limiting  Distributions  Associated  With 

Certain  Stochastic  Learning  Models’’.  Annals  of  Mathematical 
Statistics  33,  (I962). 

[3]  Karlin,  S*.  "Some  Random  Walks  Arising  in  Learning  Models  I".  Pacific 

Journal  of  Mathematics,  3>  (1953). 

[4]  Kendall,  M.G.  and  Stuart,  A;  The  Advanced  Theory  of  Statistics, 

Volume  I.  London;  C.  Griffin  and  Co.  (1958). 

[5]  Jessen,  B.  and  Wintner,  As  "Distribution  Functions  and  the  Riemann 

Z e t a  Function".  Transactions  of  the  American  Mathematical 
Society  3 8,  (1935). 

[6]  Kershner,  R.  and  Wintner,  A;  "On  a  Family  of  Symmetric  Bernoulli 

Convo lut ions . "  American  Journal  of  Mathematics ,  57 >  ( 1935 ) . 

[7]  Kemeny,  and  Snell:  Markoff  Processes  and  Learning  Theory".  Psycho- 

metrlka,  22,  (1957). 

[8]  Wintner,  A:  "On  Convergent  Poisson  Processes."  American  Journal  of 

Mathematics  57  >  (1935). 

[9]  Zidek,  J:  Asymptotic  Distributions  of  Response  Probabilities  for  a 

Stochastic  Learning  Model*  Unpublished  M.Sc.  Thesis,  Univ¬ 


ersity  of  Alberta. 


'i  .  5m.  :•  'i  ■  .0  ;  .fc  oJ  ,1  am-  '©V 


144  - 


[10]  Erdos,  P:  "On  a  Family  of  Symmetric  Bernoulli  Convolutions,"  American 

Journal  of  Mathematics,  6l,  (1939)* 

[11]  Lukacs,  E;  Characteristic  Functions,  Charles  Griffin  and  Co.  Ltd. 

(I960). 

[12]  Bush,  R.R.  and  Mosteller,  F:  "A  Stochastic  Model  With  Applications  to 

Learning."  Annals  of  Mathematical  Statistics,  24,  (1953)* 


. 


